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Abstract 

When Bose-Einstein condensates are rotated sufficiently fast, a giant vortex phase appears, that 
is the condensate becomes annular with no vortices in the bulk but a macroscopic phase circulation 
around the central hole. In a former paper [M. CORREGGI, N. Rougerie, J. Yngvason, Com- 
munications in Mathematical Physics 303, 451-308 (2011)] we have studied this phenomenon by 
minimizing the two dimensional Gross-Pitaevskii energy on the unit disc. In particular we computed 
an upper bound to the critical speed for the transition to the giant vortex phase. In this paper we 
confirm that this upper bound is optimal by proving that if the rotation speed is taken slightly below 
the threshold there are vortices in the condensate. We prove that they gather along a particular circle 
on which they are uniformly distributed. This is done by providing new upper and lower bounds to 
the GP energy. 

MSC: 35Q55,47J30,76M23. PACS: 03.75.Hh, 47.32.-y, 47.37.+q. 
Keywords: Bose-Einstein Condensates, Vortices, Gross-Pitaevskii energy. 
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Appendix B : Useful results from |CRY| 1531 

1 Introduction 

A Bose gas trapped in a magnetic potential exhibits the remarkable property that, at sufficiently low 
temperatures, a macroscopic fraction of the atoms are in the same quantum state. This phenomenon 
is referred to as Bose Einstein Condensation and has been first observed experimentally by the Jila 
and MIT groups in 1995 (2001 Nobel prize in physics attributed to Cornell, Wieman and Ketterle). 
Particularly interesting experiments consist in creating a Bose-Einstein condensate (BEC) by cooling 
atomic gases in a magnetic trap, and setting the trap into rotation. Indeed, a BEC is a superfluid and 
responds to rotation by the nucleation of quantized vortices. There is a rich literature on this subject, 
both in physics and in mathematics. We refer the reader to [A"l IFelj for extensive lists of references. 
When dealing with rotating BECs, a strong distinction has to be made according to the type of trapping 
potential that is being used. If the potential is harmonic Q there exists a critical speed at which the 
centrifugal forces overcome the trapping force and drive the gas out of the trap. By contrast, in the 
experiments of |BSSD| ISBCD| . a blue-detuned laser is applied to the condensate, resulting in a stronger 
confinement that prevents such a behavior. To illustrate the difference, let us recall briefly what is 
observed or expected when rotating a BEC at increasing speeds. 

If the trapping potential is harmonic, as in most experiments, the gas set into rotation nucleates 
more and more vortices as the rotation speed is increased. Eventually, when vortices become densely 
packed in the condensate, they form a hexagonal lattice, called Abrikosov lattice by analogy with the 
physics of type II superconductors (see [A] Chapter 5] for a mathematical study of this phenomenon). 
If the rotation is very close to the limiting speed set by the confinement it is believed that the gas will 
enter a regime where the atoms are strongly correlated. Such a phenomenon has not yet been observed, 
see }Co| for a review of the physics literature on this subject and |LSj for recent mathematical results. 

If the trapping potential increases faster than the square of the distance from the center of the trap, 
it is theoretically possible to explore regimes of arbitrary rotation rates, as first noted in Fe2J. For 
a potential of this type, variational arguments have been proposed in |FB] to support the following 
picture : three successive phase transitions should occur when increasing the rotation speed, at which 
the ground state and its vortex structure change drastically. For slow rotation speeds, there are no 
vortices and the condensate is at rest in the rotating frame. When reaching a first critical speed r2 cl , 
vortices start to appear in the condensate and organize themselves in order to minimize their repulsion. 
This leads to a vortex-lattice state where a dense hexagonal lattice of vortices is observed. There is 
then a second critical speed fl C 2 where the centrifugal force dips a hole in the center of the condensate. 
This results in a vortex-lattice-plus-hole state where the condensate is annular and still supports a dense 
lattice of vortices. When the third critical speed r2 c3 , that will be our main interest in this paper, 
is reached, vortices are expected to retreat from the annular bulk of the BEC. There should however 
remain a macroscopic phase circulation (circular superflow) in the annulus. These phenomena have 
motivated numerous theoretical and numerical studies [FJS1 IFZ1 1KB! IKF1 IKTUj . In particular, it is 
observed numerically (see e.g. |FJS| ) that in the transition regime between the vortex-lattice-plus- hole 
state and the giant vortex state, an other phase appears which displays a single vortex circle around the 
central hole. 

In this paper we aim at providing an estimate of the third critical speed in the limit where the coupling 
constant measuring the inter-atomic interactions tends to infinity (Thomas- Fermi regime). More 
precisely we build on the results of |CRYj where a rigorous upper bound to the critical speed has been 
obtained and aim at proving the corresponding lower bound. We also want to prove rigorously that 



1 which, in this context, always means 'quadratic' 
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when the rotation speed is close but below the critical speed, a circle of vortices is present in the annular 
bulk of the condensate. 

Our setting is the same as in |CDY1[ ICY1 ICRY] : we consider the Gross-Pitaevskii energy for a 
Bose-Einstein condensate in a two-dimensional 'flat' trap, i.e. a potential equal to in the unit disc B 
and +00 outside, leading to a problem posed on the unit disc. Formally this is the limit as s — > +00 of a 
homogeneous trap of the form V(r) = r s (see |CDY2| ). In the rotating frame, the energy reads 

£ GI >] := / \(V -iQx x )^\ 2 -fl 2 r 2 \^\ 2 +e~ 2 \^\ 4 , (1.1) 

JB 

and it should be minimized under the mass constraint 

|*| 2 = 1. (1.2) 



The terms in the energy have the following interpretation. The first one is the kinetic energy, including 
the contribution of Coriolis forces due to the transformation to the rotating frame. The second term 
takes into account the centrifugal force while the third one models the interactions between the atoms. 
We denote by E GP and \I> GP respectively the ground-state energy and a ground state |^] (a priori not 



unique) of ( 1.1 1 



E w := inf =£ UP [* UP ]. (1-3) 

JbI*I 2 =i 

Any GP minimizer satisfies the Euler-Lagrange equation 

- (V - iSlx^f * GP - n\H GV + 2e~ 2 |* GP | 2 * GP = / x GP tf GP , (1.4) 



where /i GP is the chemical potential associated with the mass constraint (1.2). 

To evaluate the critical speed for the appearance of the giant vortex we consider the asymptotic behavior 
of * GP and E GP as e -> 0. For |loge| < CI < e2|1 ^ ge| it has been proved |CY1 ICPRYj that many 
vortices are present and evenly distributed in the bulk of the condensate. Also the critical speed for the 
condensate to develop a hole (second critical speed in the above terminology) is estimated. It behaves as 

n ~ 2( 1 /7r£)" 1 . 

In a recent work |CRY) (see also [R] where a related problem is treated) we have proved that if 

V=^- } (1.5) 
e z \ loge| 

where CIq is a fixed constant chosen above the critical value 2(37r)~ 1 then the support of any GP minimizer 
is essentially vortex free. We have thus confirmed that a transition occurs in the regime where f2 scales 



as (1.5) from a vortex lattice plus hole state (as described in |CYp to a pure giant vortex state. We have 



also provided an upper bound to the critical speed at which the transition is expected to occur : 

f^T^ 1 (l + o(l)). (1.6) 

37T£^| logq 

In the present paper we aim at providing the corresponding lower bound. Namely, we want to show that 
if 

37T£^| log e| 



2 its existence follows by standard techniques, see e.g. |IM1I Section 2] 
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then there are individual vortices in the bulk of the condensate. 
More precisely we consider rotation speeds of the form 



37re^|loge| e^|loge| 



where f2i = fii(e) satisfies 



fii > (1.9) 
< 1 (1.10) 

The first two assumptions on f^i are natural : we want to study slightly subcritical speeds, namely we 
want to be below the threshold 2(37re 2 | loge|) _1 by a relatively small amount. The third assumption 
is technical. We do believe that the results we are going to present stay true without this assumption 
but the proofs require some improvements to apply to a regime where f2i is allowed to be extremely small. 



Stated loosely, our result is that in the regime (1.8), with £2i satisfying (1.9), (1.10) and ( |1.11[ ), there 



are vortices in the annular bulk of the condensate. They are uniformly distributed on a particular circle 
of radius i?* and we are able to estimate their number, which is of order flie^ 1 . 

If one assumes that there is a uniquely determined third critical speed separating the giant vortex phase 
from a phase containing vortices in the bulk, these results complete the proof started in |CRY) that the 
third critical speed is asymptotically equal to 2 (37T£ 2 | log e | ) 1 in the limit of small e : 

^c3 = ^-^j r (1 + 0(1)). (1.12) 

3ne 2 \ loge\ 

In fact, there exists actually a uniquely determined third critical speed. This can be seen from the fact 
that in our analysis we use the asumptions on f^i only when proving that the vortices are distributed 
along a circle. The existence of vortices in the bulk follows from our analysis under the weaker assumption 
that n ^ a(e 2 \loge\y 1 with a < 2(3+T 1 . 

On the other hand, determining the distribution of vortices for rotation speeds of the form 57 = 
a(e 2 |loge|) with a < 2(37r) _1 requires additional ingredients. As this will be the subject of a fu- 
ture work we do not elaborate more on this point here and leave aside for the present contribution the 
proof of the existence of a unique r2 C 3 . 



1.1 Statement of the Main Results 

We first recall some notation that was introduced in [CRY]. The following Thomas- Fermi energy func- 
tional will be of importance in our analysis because it gives the leading order of the GP energy when e 
becomes small with the above scaling of the rotation speed : 

f TF [p]:= / (-OVp + e-Y). (1.13) 

JB 

Here p is a matter density, it thus plays the role of | \I/ 1 2 . In particular it is positive. We minimize £ TF 
with respect to all positive p's normalized in L 1 

£ TF [ P },P>0, J b P = 1 ) 
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and find a unique ground-state given by the radial density 

p TF (r) = \ [ £ V F +eW] + = £ *f [r 2 - R 2 ] + , (1.14) 
where the chemical potential is fixed by normalizing p TF in i 1 (23) , i.e., 

/i TF = £ ,TF + £ - 2 || /9 TF||2_ (115) 

Note that the TF minimizcr is a compactly supported function, since it vanishes outside _4 TF , i.e., for 
f < Rhi where 

Rt = Jl--i T; , A TF = {R h <r <1}. (1.16) 
The corresponding ground state energy can be explicitly evaluated and is given by 

E TF = -n 2 (i ^~). (i.i7) 

We stress that the annulus A TF has a shrinking width of order e\ loge| and that the leading order term in 
the ground state energy asymptotics is — f2 2 , which is due to the convergence of p TF to a delta function 
supported at the boundary of the trap. 



When s — > with the scaling (1.51, the leading order of the energy is given by the TF energy and 
the GP density |v[/ GP | 2 is very close to p TF . In particular, the mass of <]/ GP becomes exponentially small 
in the central hole B \ A TF and the properties of the ground state and ground state energy are well 
approximated by a functional restricted to the TF annulus. For technical reasons however it is necessary 
to consider a slightly larger annulus 

A:={f : R< <r= \f\ < 1} , (1.18) 

with the choice E] 

i?< : = R h - e 8 / 7 . (1.19) 

In order to give a precise statement, we need to be more specific. Namely we need to take into account 
the leading order (i.e. E TF ) and the sub-leading order of the GP energy via a simplified functional : For 
any agZ and any real-valued wave-function cj> we introduce the 'giant vortex energy': 

£Zm--=£Ttte^-^}= f \V^ + h^+(^^-2nm-a))\^ (1.20) 



A 

where [Q] stands for the integer part of fl Let us denote g^.a and E^ p a the ground state (unique up 
to a sign that we fix by requiring that gA.a > 0) and ground state energy of this functional, see |CRY[ 
Proposition 2.3]. We recall the following result (Proposition 3.2 in [CRY]) 

Proposition 1.1 ([CRY] Properties of the optimal phase and associated density). 

There exists at least one w£Z minimizing E^ F a with respect to a. Moreover we have 



(1 + 0(| log C 1 )) (1.21) 



3 the exponent 8/7 in formula | |1.19| is choosed for definiteness : any exponent between 1 and 7/6 would do, see the 
discussion at the beginning of CRY, Section 3.1] 
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gl (n-([Q]-u;)r- 2 ) =0(1) (1.22) 

A 



where g u is the unique positive normalized minimizer of 



Note that the above proposition was proved in [CRYj with a slightly different assumption on Q, but 
the results stays true with no modification of the proof under asumptions (1.8) to ( |1.11| . 



The key to the results o f ICRYj in the super-critical case is that the asymptotics of the ground state 
and ground state energy of ( 1.1 ) are described with a very good precision by ff^ ue'™ - "^ (we will denote 



9A,u by g for short) and respectively. 



In particular, a minimizer of (1.1) contains a large vorticity similar to that generated by a giant vortex 
located at the origin, which is apparent in the phase factor e*"™ - ") 9 . To identify individual vortices in 
the annulus where the mass is concentrated we thus have to look at a reduced function as in Section 5 of 
ABM] . More precisely we extract the giant-vortex density and phase factor by introducing 

(1.23) 



■ g e i([Q]-u)e 

which is defined only on A. 

Our main results deal with the vortices of u and are stated, as is usual, in terms of vorticity measures 
of u. There actually exist two main paths to the definition of an appropriate vorticity measure. One is 
to use an 'intrinsic' vorticity measure as, e.g., in ABM, SS . This measure is defined as the curl of the 
superfluid current : 

fx := curl(m, Vu) = -curl (i (u\7u ~ uVu)) . (1-24) 

The other one, used for example in |CY1 ICPRYj . is to explicitly identify a collection of balls ^ Bj = 
B(aj, rj),j G J on the boundary of which the degree dj of u is well defined and set 

li e :=2TC^2d j S aj (1.25) 
with S a . the Dirac mass at dj. Actually, a classical technique ('Jacobian estimates' due to Jerrard and 



Soner |JS| ) allows to relate these two notions of vorticity, see Section 1.2 In this paper we will prove 
results about both the 'intrinsic vorticity measure' (Theorem |1.1[ ) and an 'explicit vorticity measure' 



(Theorem 1.2 that includes the complete definition of /u e ) : both vorticity measures are close to a Delta 



function concentrated on some particular circle of radius i?* (see Appendix A for the definition of R* ) . 

We have to deal with a technical point before stating our main theorems. In [CRY] we had identified 
the 'short-range' energetic cost (by opposition to a 'long-range' cost that we will identify in the present 
paper) of a vortex of degree 1 at some point r. It is given in terms of the following cost junction 

H(r):=\g 2 {r)\\oge\+F(r) (1.26) 



where F is defined as 



F{r):=2j g 2 (s) (cts - Qfi] - u) ^ ds. (1.27) 



4 B(a,j,rj) denotes the ball of center aj and radius rj 
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The dependence on g makes it very difficult to identify vortices close to the inner boundary of A. Indeed, 
p TF vanishes on dBf> h , so g is very small there (compared for example with its value on dB) and the 
energetic cost of a vortex too close to dBp h can not be taken into account with our method. To avoid 
this difficulty we limit ourselves to a statement about the asymptotics of /i on a smaller domain .Abuik 
where the bulk of the mass is concentrated and the density is large enough: 

Auik:={r|i?buik<r-<1}, R hn ^:=R h +e\\oge\^l\ /2 . (1.28) 
Note that, it follows from the analysis in [CRY1 Section 2.1] (the relevant result is recalled in Proposition 



B.l below) 



I* 



GP|2 



i -o{n\ /2 ) 



This justifies the notation : _4b u ik indeed contains the bulk of the mass when f2j -C 1. The radius 
R* appearing in the theorems below is an approximate minimizer of the cost function H (see Appendix A). 

Before stating the vorticity asymptotics we need to anticipate on the energy asymptotics of Theorem 
|1.3| and introduce the following notions. For any Radon measure v supported in A we define h v as the 
unique solution to the elliptic problem 



V ( iVA, 



v in A 



h v = on dA. 



(1.29) 



Note that g 2 is bounded below on A because it is the ground state of a one-dimensional Schrodinger 
operator (see Lemma 2.1 below). There is thus no difficulty in defining h v as above. Next we introduce 
an 'electrostatic energy' associated to v 



I(y) := 



A 9 



Then 



/* := inf I{u) = 1(5*) 



(1.30) 



(1.31) 



where the infimum is taken over the set T>* of positive Radon measures with support on the circle of 
radius R*. The fact that is a solution (actually, the unique solution) of the above minimization 
problem will be proved in Proposition 3.5 below. The energy (1.301 is similar to that of a charge 
distribution v in a shell with inhomogeneous conductivity described by g~ 2 . 



In the following theorems, we use the norm 



sup 

?ieCJ(.Abuik) 



U b 



V(j) 



Ja^ f' V ' 

where v is a Radon measure, to estimate jj. and /j, e . 



1/2 



(1.32) 



+ £|loge|||V0|| L ^ ( ^ bulk) 



Theorem 1.1 (Asymptotics for the intrinsic vorticity). 

Let Cli, u and /i be defined as above. Let R* be the radius defined in (A. 10), which satisfies 



R* = 




(1.33) 
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We denote S, the normalized arclength measure on the circle of radius R, 
There holds 



H(R,) . 



« — < c 

e 



H(R, 
21, 



(1.34) 



in the limit e — > 0, where H is defined as in (1.26) and I, as in (1.31). 
Remark 1.1. 



1. The vorticity measure [i is really the intrinsic vorticity quantity associated with u. Indeed, we 
expect that |u| ~ 1 (which implies curl(m, Vu) ~ 0) outside of the vortex cores, so formally, for 
any domain which boundary a does not intersect any vortex core, 



27rdeg{u, <r} ~ / (iu,d T u) 

J a 



is given by the integral of fi over the domain (using Stokes' formula). Thus, if the vortex cores are 
small, we should have [i ~ ^^2j djS a i wh ere the dj are the degrees of the vortices and Oj their 
locations. This is stated rigorously in (1.39) below. 
Note that by definition the vortices of u are identical to those of ^ GF . 



2. Let us explain a bit why (1.34) indeed captures the leading order of the vorticity. The radius R, is 
such that 

(I //;//. / " Jh 



-C- 



Taking a positive test function (or rather a sequence of test functions) in C*(.4buik) we have 



H(R, 



21, 



H(R,) 



If in addition d> is radial then 



-4buii 

H(R,) 
21, 



(r, 0)d9. 



r=R, 



H{R,) 

21* 



MR*)- 



Choosing <f> such that max0 = <p{R,) and recalling that the thickness of the annulus .4buik is of 
order e\ logej, one can obviously construct <j) such that 

<MaO>c||v^|U^ buIk)£ |io g£ |. 

If in addition the support of </> is included in a region where g 2 > C{e\ loge|) _1 we have 



1 



1/2 



< 



C £ |log £ |||V^|Uoo ( ^ bulk3 < cu\\ L ~ {Ah ^ y 



Thus 



H(R, 
21* 



which proves 



> -Ce|loge 



H(R, 



>Cn 1 |loge|||V0|U« ( ^ b111k) >C 



9.x 



1 



|V0| ; 



1/2 



(1.35) 



< C 



H(R,) 
21, ' 
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3. Let us comment on the norm ( 1.32 ) that we use in the theorem. In the denominator of the definition 



(1.321, the first term is rather natural in view of the energy estimates below. Indeed, the quantity 



ff- 2 iv</»r 



1/2 



defines a norm on C^(^lbuik) that is associated with the problem (1.31). The second term appears 
when regularizing fi in the course of the proof. 

Of course if the above norm could be controlled by e\ loge||| ^ i>\\ (A b nih) ^ or an y * es ^ function, then 



(1.341 would be equivalent to 



H(R, 
21, 



< O x | loger|, 



(1.36) 



(C c i(-4b„i k ))* 



whereas the (C^L4buik))* norm of both terms in the above right-hand side are > f2i|loge| as 
demonstrated by ( 1.35 ). However, such a control is not possible because g 2 is not uniformly bounded 
below by C(e \ loge|) _1 on .Abuik : close to the inner boundary of .4buik it is of order fi^^el loge|) _1 . 
The norm ||.|| 9 thus takes into account the effect of the inhomogeneous density background. 

4. If we define a change of coordinates mapping *4buik to a fixed annulus of radius 1 and width say 
1/2 and denote R, the image of R, under this change of coordinates and ft the push-forward of ji 
we will have 

A 



-21, 



H(R,) 



(1.37) 



in the weak sense of measures. Here 6, is the normalized arc-length measure on the circle of radius 
R,. We do not want to be too precise here in order not to mislead the reader : such a change of 
coordinates will not be used in the paper. The reason is that, applied to the energy functional, it 
would lead to a highly anisotropic (and e-dependent) kinetic energy term because of the shrinking 



width of the original annulus >4buik- The convergence ( 1.37) is useful however to identify the optimal 
number of vortices : I, is of the order of a constant, —H(R,) of the order of ilie" 1 so from (1.37) 
the number of vortices is to leading order equal to —H(R,)(2I,)^ 1 oc ilie^ 1 , as we announced. 



Theorem 1.2 (Asymptotics for the explicit vorticity). 

With the assumptions and notation of Theorem \l.l\ there exists a finite collection of disjoint balls 

B k =B{a k ,r k ),k€K 

such that Bk C »4buik and \u\ > 1 — o(l) in the limit e — > on dB k for any k G K . Defining 

fi e :=2nJ2dkS ak , (1.38) 



with 

there holds 
and 



kei< 



d k = deg(u,dB k ) 



\\^-»W(cl(A bnm )y < Cnl /2 \\oge\ « IHI ( c c i(.A bulk )) 



, H(R,) X 



« — < C 

e 



H(R,) 



21, 



(1.39) 
(1.40) 
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Remark 1.2. 



The 'vortex balls' in the above theorem satisfy some properties that we now explain. In the 
parameter regime we consider, there is a 'bad' region on which we do not have enough information 
to say something meaningfull about the vortices of u. What saves the day is that we have a control 
on the area of this region, saying that the complementary 'good' region covers the major part of 
the annulus. 

After removing fom the 'good' region a boundary layer on which the matter density is too small, 
we can cover the possible zeros of u in the remaining part by the balls in the theorem, and the area 
covered by this balls is much smaller than that of the annulus where the bulk of the mass of the 
condensate resides. 



We actually prove explicitly that most of the balls carry a positive degree and that the degree 
carried by the balls outside some region close to the circle of radius i?* is small (see Equations 



( 3.41 1, ( |3.42[ ), ( |3.43| ), ( |3.67[ ) and (44])). All this is included in the statement ( |1.40[ ). We note that 
(1.401 also implies that the vortex balls must be in some sense evenly ditributed along the circle 
Cr, , but it seems difficult to prove it without using the intrinsic vorticity measure : the electrostatic 



energy (1.30) that appears naturally in this context is not well-defined for a measure of the form 



( 1.38 ). The extra regularity of the intrinsic vorticity measure is thus used in a crucial way and there 
does not seem to be an easy way to prove ( 1.40[ ) without using (1.39), i.e. the Jacobian estimate. 
Note also that the discussion in Item 2 of Remark |1.1| implies 



thus justifying the second inequality in (1.39) 



We now present our results on the ground state energy E GF . The problem ( 1.31 ) naturally appears in 
our analysis because it is linked in a crucial way to the asymptotics of the energy E GP , as demonstrated 
in the 

Theorem 1.3 (Energy asymptotics). 

Let Qi be as above and e be small enough. There holds 



(1.41) 



in the limit s — > 0. 

The term — i?(i?„) 2 (4/*) _1 in (1.41), of order rife -2 , is a small correction compared to E^^. We 



refer to |CRYl Remark 1.4] for a discussion of the different contributions to the latter. 
The way the value — 4 f appears is through the following minimization problem 



Ah 



= inf 



|V/i„ 



A 9* 



H(R*) 



A 



(1.42) 



The above functional of v describes the energy of a vortex distribution corresponding to a vorticity 
measure v concentrated on Cn r . The right-hand side can be thought of as a renormalized energy (in 
analogy with [BBH ( SSl). The first term represents the interaction of the vortices, which is computed 
through the potential h v that they create (other formulae, including a particularly useful Green 
representation will be used in the sequel) . The second term represents the energy gain of having a vortex 
partly compensating the rotation field. The unique minimizer of (1.42) is given by — (2I :t )~ 1 H(R :t )6*, 
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which explains why the vorticity measure of u has to be close, at least in the sense of Theorem |1.1[ to 
this particular measure. Note the close analogy between our renormalized energy and that obtained in 
ABM S Theorem 1.1], the difference being essentially apparent in the weight <?~ 2 . 

We recall that superfluids or superconductors in simply-connected geometries generically exhibit 
vortex concentration around isolated points (see [A] ISS| and references therein). Here, although the 
original domain B is simply connected, the strong centrifugal forces impose an annular form to the 
condensate. As a consequence, vortex concentration along a curve occurs for this model at rotation 
rates £1 cx e _2 |loge| _1 . In [ABM], a model case for vortex concentration along a curve in the regime 
f2 oc | loge| is considered and the limiting vorticity measure is identified via T-convergence. 
Other physical situations where concentration along a curve occurs include superconductors with normal 
inclusions [AB11 IAB2j , the case described in iAAB] of a BEC whose trapping potential impose an annular 
shape in the regime f2 cx | log e | , and that of two superconductors with different physical properties in 
contact along a circle |Kac| . In those three cases the question of the distribution of the vortices along 
the curve that we solve here is left open. In particular our result is as far as we know the only one in the 
literature where the limiting vorticity measure is computed for a case where there is vortex concentration 
along a curve against an inhomogeneous background. Indeed, the inhomogeneous density profile g 2 is a 
new feature compared to the situation in ABMJ. Our method could be used in the context of AABJ 
where a similar inhomogeneity prevents from using directly the analysis of [ABM], 



A formula such as (1.12) seems to be absent from the physics literature : in |FB] the critical speed is 
estimated by comparing the energy of the giant vortex state to that of the vortex-lattice-plus-hole state. 
One then argues that the critical speed is that at which the former is smaller than the latter. This is 
too rough an estimate, for it does not take into account the fact that a circle of vortices appears at the 
transition. As a consequence, the formula we find by rigorous analysis differs from that given in |FBj 
We note that for a problem in a slightly different setting (namely, for a condensate trapped in a 
harmonic plus quartic trap) the papers |FZ( IKFj propose methods to numerically compute the critical 
speed. This is done by comparing the giant vortex energy to the energy of a condensate containing a 
single ring of vortices. Varying the parameters defining the latter state (size of the vortex cores, number 
of vortices ...) to find the best possible choice yields the critical speed. It is again defined as the one 
above which the giant vortex energy is the smaller of the two. This method does not provide an explicit 



formula such as (1.12) and of course neither does a direct numerical minimization of the GP energy. 



The response of type II superconductors to imposed external magnetic fields and that of superfluids 
to rotation of the container bear some striking similarities. This analogy between superfluidity (usually 
described by the Gross-Pitaevskii theory) and superconductivity (usually described by the Ginzburg- 
Landau theory) has been well-known to physicists for tens of years. More recently it turned out that 
mathematical tools originally developed for the GL theory could be successfully used in GP theory. 
Also there is an analogy between our (somewhat informal) terminology about critical speeds and that of 
critical fields in GL theory In particular, the analogy between the first critical speed and the field H c i 
is well-known and of great use in the papers [AABl IIMH [IM2 . 

We want to emphasize however that the Gross-Pitaevskii theory in the regime we consider largely deviates 
from the Ginzburg-Landau theory. This is due to the presence of a mass constraint (which reflects the 
fact that a BEC has no normal state to relax to, contrary to a superconductor) and of centrifugal forces 
(which could be interpreted as electric fields in GL theory). In particular, the second and third critical 
speeds we have informally defined have little to do with the second and third critical fields in GL theory 
(to our knowledge this has been first noticed in |FBJ). The discrepancy between GP theory and GL 
theory arises when the centrifugal forces can no longer be neglected, namely close to the second critical 
speed. We refer to |CDYl[ICDY2llCY] for this aspect of the theory. 



3 the latter, given in formula (20) of |FB| . is found to be rj^g —i 
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As for the third critical speed, that is our main concern here, it turns out that it bears more similarities 
with H c i than with H C 2 or H c ^. In fact, the transition happening there can be seen as a H c \ type transition 
but backwards, which is apparent in the fact that vortices disappear from the bulk of the condensate when 
the rotation speed is increased. As a consequence, many mathematical tools characteristic of the study 
of type II superconductors (strictly) between 7J cl and H C 2 will be of great use in this paper, as they have 
been in |CRY| . We refer to Section 1.3 of that paper for a heuristic explanation of this surprising fact. 



1.2 Formal Derivations 

The main intuitions in the proofs of our results are 



The rotation field being along the azimuthal vector eg with positive amplitude, a vortex of negative 
degree can only create energy. All the vortices should then be of positive degree. 



[i should necessarily have its support on Cr ( because the cost function (1.26) is minimum there. 
Indeed, the cost function gives the short range energetic cost (self-interacting kinetic energy plus 
interaction with the rotation field) of a vortex situated at a particular location. 

It is well-known that two vortices of same degree should repel each other while two vortices of 
opposite degree should attract each other. The vortices, all of positive degrees and located on Cr. 
should organize in order to minimize their repulsion. This effect leads to the uniform distribution 
on the circle. 



Let us explain how the renormalized energy (1.42 1 appears in our setting. A large part of the analysis, 



in particular the introduction of the cost function ( 1.26 1, has been carried out in [CRYj . It is reminiscent 



of the method of jAAB] although important new difficulties occur due to the different parameter regime 
we consider. 

A first step is to extract from E GP the contribution of the giant vortex profile and phase factor. This 
uses a classical method of energy decoupling originating in |LMI and an exponential decay property for 
* GP in B \ A (see [UHYl Propositions 2.2 and 3.1]). We obtain^ 

E GP = E^+£[u}+0(e°°). (1.43) 

The reduced energy functional £ is defined as follows 

4 



£[v] := / 5 2 |V«| 2 - 2g 2 B ■ (iv,Vv) + ^(1 - M 2 ) 2 (1.44) 

J A £ 



I A 

where 

B(r) := B(r)e e = (fir - ([fi] - w) r" 1 ) eg (1.45) 

and we have used the notation 

(iu,Vu) := \i(v¥v — vVv). 

The energy £[u] effectively takes into account the energy added to the giant vortex contribution E 01 ^ 
when individual vortices are present in the annular bulk of the condensate. Note in particular that in the 
case where v = 1 the energy is 0. We are thus interested in estimates for £[u]. A strictly negative value 
of this energy will indicate the presence of vortices. We note that this energy functional is very similar to 
a functional appearing in |AAB| where an annular condensate at slow rotation speeds is considered (see 
also |IM1[ IIM2| ). The major difference is that the domain A depends on e in a crucial way : its width 
tends to zero proportionally to e\ loge] when s — > 0. We refer to |CRY1 Section 1.3] for a more detailed 
discussion on this point. 



6 The notation 0(e°°) refers to a quantity decreasing faster than any power of e, e.g. exponentially fast. 
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The study of the energy £[u] starts with the construction of vortex balls : we isolate the set where u could 
possibly vanish and cover it with a collection of disjoint balls with relatively small radii. The growth 
and merging method, introduced by Jerrard and Sandier independently (see |SS[ Chapter 4]) yields lower 
bounds of the form 

' ff 2 |V«| 2 > Trldilfl^Ca,)! loge| (1.46) 

Bi 

where B t = £?j(ai,7"j) is a generic ball in the collection and di the degree of u around the ball B^. We 
neglect remainder terms for the sake of clarity. We stress that strictly speaking we do not cover all 
possible zeros of u by this method : there is first a layer close to the inner boundary of A where the 
density g 2 is too small to use this method. More importantly, the shrinking width of A makes it difficult 
to obtain the estimates needed for the implementation of the method. There is thus a region where we 
have basically no information on the vortices of u. We will neglect this fact in this sketch since it produces 
only remainder terms in the energy. 

The next step is to use a Jacobian Estimate (see |SS[ Chapter 6]) to obtain 



(1.47) 



where 5 a . is the Dirac mass at a,. We then note that the function F in (1.27) is constructed to satisfy 



V ± (F) = 2g 2 B, F(i?<) = 0. 



Thus, integrating by parts the second term in (1.44) we have 



2g i B-(iv,Vv)= Ffj,- F(iu,d T u) 

A J A JdB 



(1. 



The boundary term has been shown to be negligible in CRY . At this stage we thus have essentially, 
gathering Ol, fl48| and |L47|, 



/ g 2 \Vu\ 2 -2g 2 B-(iu,Vu)>Y^2n(l\d l \g 2 (a l )\\oge\+d l F(a l )) + f g 2 \Vu 

J A i \ z J JAXUiBi 



(1.49) 



In |CRY| the first two terms in the above right-hand side were enough to conclude. Indeed, F is negative 
in the bulk, thus we can already read on the above formula that negative degree vortices can only increase 
the energy. All vortices should then be of positive degree and we can write 



/ g 2 \Vu\ 2 -2g 2 B ■ (iu,Vu) > V 27rd t ff (a 4 ) + / g 2 \\7u 
J A i JA\U,Bi 



(1.50) 



where H is the (radial) cost function (1.26). Above the critical speed, H is positive everywhere in the 



bulk, thus no vortices can be present in the condensate. By contrast, when the rotation speed is sub- 
critical, the cost function has a negative minimum at i?* (actually, close to R*, see Appendix A for the 
precise definition of R*). Positive degree vortices are thus favorable in this regime, and they will gather 
close to the circle Cji t . 

We have now to understand what mechanism limits the nucleation of vortices and leads to their uniform 



distribution on Cr„. This is where the third term in the right-hand side (1.49) comes into play. In [CRY 
this term was neglected, but in the present context it is crucial to bound it from below. 
Let us introduce the superfluid current of u 



j := (iu,Vu). 



(1.51) 
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Note that by definition fj, = curl(j'). The idea is now that the kinetic (first) term in (1.44 1 is essentially 
due to the presence of this supercurrent. The energetic contribution of the current inside vortex balls is 
taken into account in (1.46). The next step is to estimate the energetic contribution of the part of the 
current that is located outside vortex balls (this is the long-range energetic cost of the vortices we were 
alluding to before). To this end we define a modified superfluid current 



Jmod :_ 



jm.A\ UiBi 
in U, Ba, 



and the associated vorticity 



Mmod := curl(j mod ). 
It turns out that (in a sense that we keep vague in this sketch) 



Mb 



Then, combining (1.47) and (1.54) we can rewrite (1.50) in the form 



g 2 \Vu\ 2 - 2g 2 B ■ (iu,Vu) > / H/j, 



A\UiB t 



g 2 \Vu\ 



(1.52) 

(1.53) 
(1.54) 

(1.55) 



Let us now estimate the last term in the above expression. By definition of the vortex balls, outside the 
set UiBi we have lul ~ 1, thus 



7 2 |Vd 2 



-4\UiBi 



Recalling the equation (1.29) satisfied by 

curl I j mod ■ 



A\UiBi 

we have 



9 2 \J\ 2 



= in A. 



(1.56) 



If we were working in a simply connected domain this would imply that there exists a function / such 
that 

Jmod = -4V X /l^od + V/. (1.57) 



The argument has to be modified because we work on an annulus, but since this does not produce 

14, we 

(1.58) 



significant modifications of the energy, we assume that (1.57) holds. Then, using /i Mmod = on dA, we 
obtain 

' ' ' ' |2 



f 5 2 |V W | 2 ~ f 5 2 |j mo d| 2 > / 4>I V/ V„„ d 

JA\UiBi J A J A 9 



Gathering (1.55) and (1.58) we obtain, up to the (many !) remainder terms we have neglected 

1 



g 2 \Vu\ 2 - 2g 2 B ■ (iu,Vu) > 



A 9' 



\Vh, 



(1.59) 



The right-hand side of the above equation is equal to the rcnormalized energy of /i mo d defined in (1.42) 
if this measure has its support on Cr„- Intuitively this is justified because the vortices of u will want to 
be close to Cn t where the cost function is minimum. Mathematically, showing that the right-hand side 
of (1.55) is bounded below by the infimum in (1.42) uses essentially two ingredients. One is the rigorous 
version of the Jacobian Estimate making precise the property (1.47), the other is the information that 
the potential associated (via (1.29)) to the measure minimizing (1.31) is constant on the support of the 
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measure. This is a classical fact in potential theory. In our context this means that hg, — ft* is constant 
on Cr,- We refer to Section 3.3 for details on the way we use this information to complete the lower 
bound, and accept for the present sketch that /i mo d has support on Cr, in the limit e — > so that (1.591 
yields the renormalized energy (1.42) of ^ mo d as a lower bound to £[u\. It is then intuitively clear, if 
we believe not to have lost too much information on the way, that /i mo( j should 'almost' minimize the 
renormalized energy. Thus it should be close to — i?(i?*)(2/*) _1 <5* and by (1.54) so should \i. 



What remains to be proved is that 'we have not lost too much information on the way', i.e. we want 
an upper bound to the energy matching our lower bound. We thus construct a test function displaying 
the optimal number of vortices evenly distributed on the circle Cr„ . To evaluate the energy of such a 
test function we adapt the well-established technique used e.g in jAABi lASSi lABll IAB21 IABM[ [55]. 
It involves in particular a Green representation of the electrostatic energy (1.30). In our case the 
computation is significantly complicated by the fact that the matter density g 2 vanishes at the inner 
boundary of A and more importantly by the particular geometry of the annulus (fixed radius and 
shrinking width). We have to rely on the periodicity (in the angular variable) of our test function to 
obtain the estimates allowing to conclude that the appropriate upper bound holds true. 

The paper is organized as follows. We first construct our trial function and prove our upper bound 
in Section [2] Section [3] is then devoted to the proof of the lower bound. We will collect many estimates 
on the way that will allow us to conclude the proof of Theorems |1.1| and L2 in Section [3] Appendix A 
is concerned with the analysis of the cost function and contains mainly technical computations whereas 
Appendix B gathers several useful results from [CRYj that we recall for the convenience of the reader. 



2 Energy Upper Bound 

In this section we provide the upper bound part of the energy asymptotics of Theorem |1.3| by constructing 
a trial function for £ GP which displays the expected vortex structure. 
The result we are aiming at is the following : 

Proposition 2.1 (Upper bound to the energy). 

Let Af = Af(s) be a number satisfying 



(2.1) 



in the limit e — > 0, where 2* is defined in (1.31) and H is the cost function (1.26). We have, for e small 
enough 



£ GF < E% p u + 47r 2 (A0 2 (l + o(l))/, + 2-kNH{R*) + O 



/ AHogllogel 
V e|loge| 



E 



gp H(R*) 2 



47, 



(l + o(l)). (2.2) 



2.1 The trial function 

We start by dividing A into N oc j ^ ge | identical cells. The cell Ai is defined as 

A l =AC\{8 l <e<e i+l \ (2.3) 

with 

9i = (i-l)—, i = l...N. 
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Figure 1: Vortex configuration in the trial function. 



and the convention that 9n+i = We have that 9i+i — Qi is independent from i (that is the cells are 
identical) and proportional to e \ \oge\ (which is the thickness of the annulus ^4). We construct a structure 
of vortices periodic in the angular variable, of period 8 i+ i — 6i = jf. 

We introduce a parameter t which will be the size of the vortex cores and M points . • . ,Pi.M in 
each cell, which will be the locations of the vortices. We require that for any i, j 



B(p id ,t) C Ai 



(2.4) 



and that the collection of points {j>ij)i=i...N, j=i...M is evenly distributed on the circle of radius i?* (which 
we denote CrJ. The distance between two adjacent vortices is thus 2i?* sin(7r/MiV) ~ 2nR*(MN)~ l . 
The number of vortices will be fixed (i* is of the order of a constant, see (3.55), H(R^) is negative and 
proportional to flie^ 1 , see (A. 13) and ( |A.18 )) so as to satisfy, in the limit e — > 0, 



N := MN 



471-7* 



(l + o(l)) ex 



fin 



(2.5) 



The natural vorticity measure associated with our ring of (degree one) vortices would be a sum of Dirac 
masses 2n ^ S Pi . , which would lead to an infinite energy. We thus introduce a regularized vorticity 
measure 



^ Id B(p«, J ,t) 



(2.6) 



where Ids( Pi ,t) is the characteristic function of the ball B(pij,t). The measure / is normalized to have 
total mass 2ttMN (i.e. 2n times the number of vortices). We require in our construction that this 
function be invariant by some reflections of the angular variable. More precisely, for any i, we require 



f(r,6) 
f(r,6) 



f (r, 26i 
f{rA- 



(2.7) 
(2.8) 
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Clearly this conditions imply the periodicity of / : 

Mr, 0) := f\ Ai (r, 6) = f ]Aj (r, 6 + 6,- 6,). (2.9) 

Note that the freedom that we have in the choice of the number of vortices (in the o(l) term of formula 
(??)) is used to ensure that it is possible to construct a vortex configuration having the symmetries we 
want. A typical distribution of balls guaranteeing these symmetries is one that is even (in the angular 
variable) in each cell with the same distances of the two extremal balls to the radial boundaries of the cell. 



Our trial function is given by 



* = I 96 VmA (2.10) 

1 in B \ A 



with 

u = c£e^. (2.11) 
The real function £ is a cut-off function ensuring that 

• v vanishes in small discs (of radius t) around the vortex locations, where the phase <fi will have a 
singularity 

• v vanishes continuously at the inner boundary of A, so that ^ is indeed in the energy space. 

in such a 



The constant c is chosen so that J B \^\ 2 = 1. As for the phase (j>, we define it in Lemma 
way that v has degree one around each vortex in the collection. 



2.2 



Let us introduce functions £j j ,i = l...N,j = l...M satisfying 

1 hxA\B( PiJ ,2t) 



ai-l?! 117 !^ ^ (2-12) 



and 

Obviously one has 
and we can impose 



< i,., < 1. 

V&j = in B{pij,t) U B{ PiJ ,2t) c (2.13) 
iV&il < j in B(p itj ,2t) \ B(pij,t). (2.14) 



Next we define, for some radius R larger than 

f 1 if r > R 

^ = \^Tifr<R. (2 ' 15) 

I, R—R < 

The subscript " BL" stands for boundary layer and is justified because our choice of R will satisfy 

R = i?< + e n < R h - ^e 8/7 (2.16) 

for some large enough power n (we will use later the fact that one can choose n as large as one pleases 
but for the present n — 2 is sufficient) and s small enough. The inequality above is a consequence of 
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( 1.19 1 and allows one to use the exponential smallness of g 2 (proved in [CRY] and recalled in Proposition 
B.3I in the region i?< < r < R. 
The final cut-off is given as 

£ iiu.HiV (2.17) 

The choice of £bl vanishing continuously at r = i?< will ensure that \& £ H 1 (B). We now choose 

t = e^ 2 \\oge\ 1 / 2 . (2.18) 

This is the optimal choice and allows to estimate the L 2 norm of ^ : it is easily seen that for any i and 
any j 

g 2 = 0(e 2 ), 

B(Pi,j,t) 



because g 2 < C(e|loge|) . Moreover, using the exponential smallness result (B.4| for g 2 one can see 
thafl 

/ g 2 = o(e°°). 

jR<<r<R 

With these two estimates in hand, using the normalization of g 2 we conclude that 

g 2 £ 2 > 1 - CMNe 2 . (2.19) 



On the other hand 



9 2 e < g 2 = i 



I A J A 

and thus one can normalize "J/ by taking a constant c satisfying 

c 2 = l + 0(MNe 2 ). (2.20) 

We now turn to the definition of the phase of v. The following modified density p, defined for r G A, 
will be used in the definition: 

V F (r) if r > R 
g 2 (r) if r < R 



where 

R < <R = R h + e 5/6 



with the actual choice specified in Lemma 2.1 below. This function is constructed to satisfy two properties 
that we shall need in the sequel. First its gradient has to be bounded by C(e\ loge|) -2 for r > R, which 
is insured by the explicit form of p TF . This is the main reason why we do not use simply g 2 : the bound 
that is available on the gradient of this function CRY, Proposition 2.7] is not sufficient for our purpose. 
Second p has to stay close to g 2 in L°° norm, which is insured by the following 

Lemma 2.1 (Properties of p). 

Let us define 

R:=R h +e 5/6 . (2.22) 

Then for any r G A 

\g 2 (r)-p(r)\<Ce 3 / 4 \loge\ 2 p(r). (2.23) 
Moreover p is bounded below by a positive constant in A (the constant depends on e). 



7 saying that some quantity is 0(e°°) means that it goes to zero faster than any power of e when e — > 
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Proof. The estimate is of course t rivia l when r < R. When r > R it is a consequence of Proposition 2.6 
in [CRY], recalled in Proposition B.4 and of the explicit formula for p TF . The fact that p is bounded 



below is a consequence of the corresponding result for g . Indeed, g satisfies the equation 



-Ag 



([ft] 



2„3 



-g-2Cl([Q] -u)g + 2e- 2 g- 



p gp 9, 



(2.24) 



with Neumann boundary conditions on OA, where p, is the Lagrange multiplier associated with the 
mass constraint. It is well known that such a function, ground-state of a one-dimensional Schrodinger 
operator, cannot vanish except at the origin (see e.g. [LL( Theorem 11.8]). A possible proof is through 
the Harnack inequality [Moj . □ 



Note that p has a jump discontinuity. We could have constructed a regular function instead, but since 
the discontinuity has no consequence in the sequel, we stick to the simple definition (2.21 ). 
We now define hf as the unique solution of 



V ( iV/i 



hf = on OA. 



fin A 



(2.25) 



Note that the elliptic operators app earin g in (1.291 and (2.25) arc similar up to the replacement of the 
weight g~ 2 by p. Thanks to Lemma 2.1 these two operators are close in some sense. We denote 



ldhf 
,,B P dn 



2tt 



1 

2^ 



ldhf 
,)B P dn 



(2.26) 



where #; denotes the deritaive in the outward normal direction. Also we introduce T as the solution of 



-V ( iVL ) = in A 



dn 



1 on dB 
on dB R< . 

Note that T is radial and has the explicit expression 



(2.27) 



(2.28) 



The denominator in the above equation is a 0(1) thanks to (2.23) and the normalization of g 2 . 
We prove the following lemma, which defines the phase of our trial function. 



Lemma 2.2 (Phase of the trial function) 

The formula 



= -V x h f 
P 



where 



ht = ht — 



defines a phase, i.e. e 1 ^ is well- defined, in A\UijB(pij,t). Moreover 



PC 2 |V0| 2 < 



1 

A P 



Vhi 



< 



Nh, 



A P 



•0(1). 



(2.29) 
(2.30) 

(2.31) 
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Proof. By definition of hf and T, the right-hand side of (2.291 is irrotational in A \ UijB(pij,t), thus </> 
is well-defined locally in this set. To see that e 1 ^ is well defined, we must check that for any closed curve 
o~ included in A \ Ui,jB(pi,j, t) 

V0-tg27tZ. (2.32) 



Clearly, it is sufficient to consider two cases : 

1. (j winds around at most one ball B(pij,t) 

2. cr is a contour enclosing all the balls 

and one deduces the general case from these two. 

In case 1, (2.32) is a simple consequence of the quantification of the mass of /. One shows easily that 
(2.32) is satisfied using (2.25), (2.27) and integrations by parts. In case 2, equations (2.25) and (2.27) 
give 



• T 



ldh. 



1 dT 



p dn I |vr| 2 P dn 



But, using again (2.27) 



1 dT 

,,B P dn 



,>B P dn 



- ivri 2 

A p 



and thus it is clear from (2.29) that 



ldh j 



1 dT 



OB 



p dn J A 1 | vr| pdn ^ 



€ 2?rZ. 



To prove (2.31 ) one remarks that 







2 






Vh f 




f l 


Ja p 




Ja P 



Vh 



f 



Ia 1 p \^ 



rvr 



Then, using the equation for T and the boundary condition on hf 

1. 



A P 



Vhf • vr = 0. 



Finally one uses that < n < 2tt and 



- |vr| 2 > c 

A P 



which follows easily from (2.28). 



We can now start the computation of the energy of our trial function : 



(2.33) 



□ 



Proof of Proposition \2. 1\ A classical technique of energy decoupling (see e.g. the proof of Proposition 
3.1 in CRY]) uses the variational equation satisfied by g to show that 



■S[v] 



where 



E[v]:= / g 2 \Vv\ 2 -2g 2 B-(tv,Vv) + ^(l-\v\ 2 ) 2 



(2.34) 
(2.35) 
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and 



B(r) := B(r)e B = (Qr - ([fi] - u>) r" 1 ) 



(2.36) 



We have also used the notation 

(iv, Vv) := \i(vWv — vVv). 

We now have to estimate the energy £[v}. We first note that it can be decomposed into the contribution 
of the density profile £ and the contribution of the phase, i.e. the energy generated by the vortices : 



£\v] = c 2 



A 



(l-c 2 C 2 ) 2 )+ ^(cV^ 2 |V^| 2 -2 ff 2 S.(z«,V W )). (2.37) 



The contribution of the profile (first term in the above equation) is readily estimated by using the 
definition of £. We separate the boundary layer {r < R} where £ = £bl from the bulk where £ = J\ i ^ £jj 



5 2 iv£i 2 + 9 - 2 (i - c 2 e 2 ) 2 ) = / Uwe + 4 (i - c^ef 

A \ £ J JAn{r<R} \ £ 



+ 



An{r>R} 



<7 2 |V£| 2 + ^(l-c 2 £ 2 ) 2 



In the boundary layer we use (B.4) and the definition of £bl (2.15) to obtain 



3W + ^(l-c 2 £ 2 ) 2 ) =0(0- 

An{r<R} V 



For the bulk term we observe that, since g 2 < C(s\ loge|) , (2.13) and (2.14) imply 

MN 



I ff 2 |V£| 2 <C 

JAn{r>R} 



e|loge| 



We then write 



/ 



>R} 



1-c^ 



An{r>R} 



, I 9V-~ f 9 4 e) (2.38) 

JAn{r>R} £ JAn{r>R} J 



We remark that £ = 1 outside of UijB(p iy j,2t), thus 



/ . 

JAn{r>R} 



e|loge| 



using (2.18) and g 2 < C(e\ loge|) 1 . For the second term in the right-hand side of (2.38) it is sufficient 



to recall that g 2 < C(e|loge|) 1 , £ < 1, \A\ oc e|loge| and use (2.20) to obtain 



fc 2 -l) ( , 



MN 



£ -IAn{r>R} £ JAn{r>R} ) £ \lOg£\ 



We can thus conclude that (using again (2.20)) 



M < (1 + 0(MNe 2 )) f ( 5 2 e 2 |V^| 2 - 2g 2 B ■ (iv, V«)) 



+ C 



MN 
e|loge| 



(2.39) 
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We thus need to compute 



A 



which is the term taking into account the kinetic energy of the vortices, and 

-2 / g 2 B ■ (iv,Vv) 



(2.40) 



(2.41) 



which represents the interaction of the vortices with the rotation field. These two terms will be estimated 
in the next two subsections respectively. To conclude the proof it is sufficient to gather (2.341, (2.39), 
(2.42) and the results of Proposition 2.2 and 2.3 below. 

□ 



2.2 Kinetic energy of the vortices 



We first relate (2.40) to the energy of hf. using (2.23) and (2.31) we have 

|. 9 2 e 2 iv^i 2 = (i+o( e 3 / 4 iio g ei 2 )) J P em 2 



< fl + 0(e 3/4 |loge| 2 )) / - |V^-| 2 + O(l). (2.42) 
v ' JaP 



The required estimate is then given in the following proposition 

Proposition 2.2 (Kinetic energy generated by the vortices). 

Let Af be as in (2.5) and I* defined by (1.31). For e small enough there holds 

( A/" log | log e| 



A P 



- |V/i/| < 4^ 2 (A0 2 (1 + 0(1))/* + irNg 2 {R*) |loge| + O 



V e|loge| 



(2.43) 



Note that the third term in the right-hand side of the above equation is always much smaller than 
the second one because g 2 (R*) oc (e| loge|) _1 . It is also a remainder compared to the first one because 
of the specific choice (2.5) and the assumption ( |1.11[ ). 

The proof requires several constructions. The term we want to compute can be estimated using a Green 
function defined as follows 



(2.44) 



{^x) V ^ G ( x ^y)j = S v( x ) for xeA 
G(x,y) = for x £ dA. 



The existence and symmetry (G(x,y) = G( y, x) ) of such a function is classical, because at fixed e, p is 
bounded above and below in A (see Lemma 2.1). Then, using the Green representation of hf 



1 

A P 



G(x,y)f(x)dxf(y)dy. 



(2.45) 



IxGA JyeA 

Before going to the technical implementation, let us explain the ideas behind the computation. We have 
/ / G{x,y)f{x)dxf{y)dy=Y J II G(x,y)^-. (2.46) 

JxeAJyeA i.j,k,l J J B(p ii3 ,t)y.B(p k , L ,t) 1 

As is well-known, the Green function has a logarithmic singularity at x = y. More precisely, for x close 
to y we expect the behavior (see Lemma 2.4 where this is proved for slightly different Green functions) 

g 2 (y) 



G(*, y )~^bg 



ejlogej 
\x - y\ 



2tt 



log 



el log el 

\x - y\ 
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The second estimate is a consequence of Lcmma pO] We replace G by the above expression in the diagonal 
terms (i.e. i = j and k = I) of the sum (2.46). A direct computation yields the second term in (2.43). 



For the off-diagonal terms we simply note that G is regular when x is not too close to y, in particular 
when x G B(pi,j,t) and y € B(pkj,t) with i ^ k or j ^ I. Then the off-diagonal sum can be seen as a 
Riemann sum and yields to leading order 



£ // G ^y) 



Adxdy 



(j,i)/(fc,0 



{AttNY I G{x,y)5*{dx)5*(dy) = (AirAf) 2 -\Vh s J 
>AxA ja P 



There remains to note that because p ~ g (Cf Lemma 2.1 ) 



A P 



A 9 



:|V/i 4 . 



I*. 



For this heuristic presentation we have deviated from the actual proof procedure that we present below. 
We note that sever al aut hors (see e.g. |AAB1 IASS| IAB11 rA"B2l IABM( IS5 ] ) have already successfully used 
the representation (2.45) for the computation of similar quantities. However, in our case, the particular 



geometry of A (fixed radius but shrinking width) makes it difficult to obtain the properties of G required 
in the computation. We thus prefer to take advantage of the symmetry of the vorticity measure / to 



obtain another expression of (2.45) 



Let us introduce for any i the function hj defined on the cell Ai as follows 



'-V(jV/i}) =fi in Ai 
K) = on dA l n dA 
i^f =0on 8Ai\dA. 

p on 1 \ 



(2.47) 



See (2.9) for the definition of /j. Note that we impose Dirichlet conditions only on the azimuthal part of 



the boundary of Ai (that is, the part that coincides with the boundary of A), while the behavior on the 



radial parts of the boundary is left free, which leads to the Neumann condition in (2.47). Because of the 
symmetry of /, the following holds 



Lemma 2.3 (Alternate definition of hf) 

We have 



h) on A 



(2.48) 



i.e. for any i £ {1, . . . -ZV} 
As a consequence 



hf = h % f on A 



-iw 2 

A P 



f 



(2.49) 



Proof. Using ( 2.8 1 and the fact that p is radial, we obtain that hj is symmetric with respect to a reflection 
of the angular variable, that is 

h f (r,6) = h f (r,6i + 6 i+1 -6). (2.50) 



On the other hand, using (2.9) we have also 

h)(r,6) = h f (r,8 + l 



(2.51) 
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Using the last two properties, it is clear that J^i ^} is well-defined, continuous on A and satisfies 



< +00. 



A simple computation shows that 



P 




on A, in the sense of distributions. The conclusion follows by uniqueness of the solution to the elliptic 
problem ( |2.25| ). □ 

The point of using this expression of h f is that we are now reduced to the computation of the terms 
in the right-hand side of (2.49|. We introduce new Green functions, defined on each cell 



S y (x) for x £ At 



Gi(x, y) = for x £ dAi D dA 
= for x £ dA l \ dA. 



(2.52) 



1 dGj(x,y) _ 



Then, using the Green representation of h\ 

1 



A, P 



m\ = 



Gi(x,y)f(x)dxf(y)dy. 



(2.53) 



Now we are reduced to a computation on each cell. What makes the computation easier compared to 
what it would have been using (2.45) is that each cell can be mapped by a dilation of the variables to a 
fixed domain. The difficulty of the shrinking width of A is thus avoided. 

Note also that instead of computing in (2.45) the interaction through the Green potential G of each 
pair of vortices in our collection, we can now simply compute the interaction through Gi of each pair of 
vortices in Ai, and sum the contributions of each cell. Thus each vortex effectively interacts only with 
M oc fi]J loge| vortices through a modified Green function instead of NM oc fiie -1 vortices through the 
original Green function. 

We now prove useful properties of the Green functions defined in (2.52 ). It is mostly here that our choice 
of p will prove useful, in particular the fact that the gradient of this function is properly bounded above 
close to Cfi t . Since / has its support close to the circle of radius i?* we are mainly concerned with 
properties of Gi in that region for the computation of (2.53). 

Lemma 2.4 (Properties of the Green functions). 

Let Gi be the Green function defined in (2.52). Let Cq be some constant, small enough and independent 
of e. There holds 

1. for any y G A,, G(. , y) € W 1 ' p (A i ) for any 1 < p < 2 
2- Gi(x, y) > for any (x, y) £ A t x A. t \ {x = y} 
3. Gi is symmetric, i.e. Gi(x,y) = Gi(y,x) 

y\ — R*\ < C(3£|log£|. For any compact set K CC Ai whose diameter 



Let y £ Ai be such that 
satisfies 

diam(iT) < C G e|loge|, 
and any y £ K , there exists a constant Ck such that, for any x £ K 



Gi(x,y) 



p{y) 

2tt 



log 



\x-y\ 

e|loge| 



< C 



K 



1 

e| logs) " 



(2.54) 



(2.55) 
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5. for any y G Ai, Gi{x,y) is continuous, uniformly in e, on 

An{|H-i^| <C G e|loge|}n{|a;-y|>C G e|loge|}. 

Before giving the proof of this result we stress that we are going to use some results of [St] (see also 
LSWJ), where the main emphasis is on homogeneous Dirichlet boundary conditions. However, a careful 
inspection of the proofs shows that the results apply as well to our case where the boundary conditions 
are mixed (Neumann/Dirichlet), as long as we are only concerned with interior estimates, as we are. 

Proof. The first three properties are classical, because the operator we consider is uniformly elliptic in 
Ai (note however that the ellipticity is not uniform with respect to e). 
To prove items 4 and 5, we introduce a dilation of the variables 



mapping Ai to (e| log; 



x = e\\oge\x, y = e\\oge\y 
~ Y Ai, which is a domain of fixed size. We define 



p(x) 
G(x,y) 

w(x,y) 



e\\oge\p(x) 
e\\oge\Gi(x,y) 

G(x,y) + ^log\x-y\ 



We have from (2.52) 



and 



-V s ( j^VzG(x,y) ) =S s (x) 



— — V £ w(i,y) 
p{x) 



V £ (log|x-y|) :=F(x,y) 



(2.56) 

(2.57) 
(2.58) 

(2.59) 

(2.60) 
(2.61) 



2irp 2 (x) 

on (e| logel)" 1 ^ in the sense of distributions. We now proceed as in [ASSl Lemma 3.1] : 

Let y £ Ai be such that \\y\ — i£*| < C G e| loge| for a small enough constant Cq and if be a compact set 

included in Ai satisfying (2.54) and y G K. We have, for any x e K 

\\x\-R,\ <2C G£ |log £ | (2.62) 



for a small enough constant Cq. Let K be the image of K under the dilation of variables (2.56). Using 



Theorem 2 of [Me , there exists p > 2, p' < 2 and a constant Ck such that (here we fix y and consider 
the functions as depending only on x) 



\ LP( K) < C'k (\\V w \\ L p'(K>) + ll^llw-i.: 



(K'\ 



(2.63) 



where K' is some set such that K C K' . Using a Sobolev embedding, we can take some 1 < q < 2 such 
that 



l^llvr-^CK-') - c \\ f \\li(k>) 



(2.64) 



But, using the definition of p and the explicit expression of p TF we see that as long as x satisfies (2.62) 
one has p{x) > C(e\ loge|) _1 and thus 

p(y) 



p 2 (x) 



< C 



because p(y) < C(e\ loge|) 1 for any y e A. Also, if Cq is small enough , p(x) = p 
and thus 

|Vp(x)| < — ,2 

e z \ logej^ 



TF , 



x) for any x £ K , 
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on K, which implies that on the domain we are considering 



2np 2 {x) 



< C. 



Hence, for any q < 2 

\\F\\ Lq{Kr) < C. (2.65) 

On the other hand, using Theorems 9.1 of [St] (note that the value of the constant in the right-hand side 
of (9.6) in [St] is given by Theorem 4.2 of the same paper) to estimate the W 1,p norm of G we have, for 
any p' < 2 



\\Vw\\lp\k>) ^ \\ w \\ w i.p'(K>) < C | (e| logej) 1 A I 



,1/2-1/p' 



(e|lo ge |)-i^ 



Sy < C 



(2.66) 



because (e|loge|) l A% has a fixed size. Plugging (2.64), (2.65) and (2.66) in (2.63) and using a Sobolev 
embedding (recall that p > 2) we obtain 



IHIl<~(AT') < C\\w\\ W ^p{K') ^ C - 



(2.67) 



There only remains to change variables to conclude that (2.55) holds. 

To prove Item 5 it is sufficient to fix y £ (e| loge|) _1 ^li and show that G(x,y) is continuous on 

(e| \oge\)~ 1 A l n - (e\ logeD" 1 ^ < C G ) n {\x - y\ > C G } . 

This follows from Theorem 7.1 of [St]. Note that the constant in the right-hand side of (7.5) in that 
paper does not depend on the domain, so that the continuity is indeed uniform with respect to e. □ 



We are now able to present the proof of Proposition |2.2[ With the above constructions in hand there 
is now no real difficulty to adapt a method that has already been used several times in the literature 
[£SSl rXBll IXB21 KBMlISS] . We will thus be a little sketchy in places. 



Proof of Proposition 2.2. Step 1. Clearly, from ( 2.51 1 we have for any i,j £ {1, . . . , N} 

Ja % P J Aj P 



It is thus sufficient to compute one of the terms in the sum (2.49) and the result will follow. 



We use the Green representation of for the computation : 

Gi(x,y)f(x)dxf(y)dy 



-\m\ 2 = 



Ai P 



(x,y)£AixAi 



E 



(i,!()eB(pi 1 i 1 {)xB(j)i, j ,i) 



Gi(x,y) — dxdy 



+ 



E 



i#i,|p l>4 -p l>J |< j 8 e |l ge|" / " / (*'») eB (*'M,*)xB(pi,i,t) 

Gi(x,y)f(x)dxf(y)dy 

(ar,y)e^iX^in{|a;-i/|>^e|loge|} 



G\{x,y)^dxdy 



(2.68) 



where /3 is a small parameter (see below). On B(pi i, t) x B{p\^ 1 1), according to (2.55 ) there is a constant 
(independent of i) so that (recall that t — e 3 / 2 | logej 1 / 2 and \p\^\ = R*) 



\x-y\ 
e|loge| 



< 



C K 
e|loge| 
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Also, on B(pi,i,t), p = p TF which gradient is bounded by C(e\ loge|) 2 , thus on B{p\,ii t) x B{pi^, t) 

C K 



Gi{x, v) + — — log — 

2ir \e loge 



< 



which allows to compute 

E 



Using ( |2.18| and ( |2.23[ ) we thus obtain 



Gi{x,y)-^dxdy = 2nMp(R* ) log 



e|loge| 
e|loge| 



ill 



E 

i 

We recall that 



(£c,y)eB(pi, i ,t)xB(pi ii ,t) 



G 1 (x,y)^dxdy = nMg 2 {R*) \loge\ + O 



e|loge 

/ M log | log e 
V £•! log e| 



(2.69) 



(2.70) 



|pi,i -Pl,i+l| = 2i?* sin(7r/MiV) - > t. 



Combining this fact with the estimate (2.551, the upper bound p < C(e|log£|) 1 and a Riemann sum 
yields 



E 



&i,\Pi*-Pl4\<Pe\lKer J (*' y)eB te.*> t )x B fa.i> t ) 

c 



Gi(x,y)- i dxdy 



< 



e\ loge 



E E 



» j^h |Pi,4— Pi,il<i8e| loge| 



lof 



e|loge| 



log 



|»-P<| 



where we denote 

Cjj, := {leCfl.^ifl^MiV)- 1 < lar-ftl </3e|loge|} 
We conclude from ( |2.7l| , using MiV cx ili£ _1 that 



e|loge| 



1 (2.71) 



E 

j#«> |pi,i— Pi,j I </3<r| loge | 



(i,y)eB(pi,i,t)xB(pi,j,t) 



Gi (*, y) -dxdy < Cj9| log 0| A/ 2 TV. (2.72) 



The term on the fourth line of (2.68) is also estimated using a Riemann sum. Using the fact that G\ is 
continuous uniformly in e on Ai x Ai\{x — y} we obtain 



II(x 



(x,y)eAixAin{\x-y\>/3e\ loge|} 



Gi(x,y)f(x)dxf(y)dy 



4tt 2 (MiV) z 



(x,y)£A 1 xA 1 r\{\x-y\>0e\ loge|} 



d^y)^^^) (l + o(l)) 



< 4tt 2 (MNY 



(x,y)eAixAi 



G 1 (x,y)6 if {x)d if (y)(l + o(l)). (2.73) 



Let us denote h\, the function defined in Ai satisfying 



■V 



= (5* in A 



h\ = on &4 4 n dA 

ldhl 
. P On 



(2.74) 



=0on &4i\d,4 
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and ft,* satisfying in A 

f-V(iVJi*) =S« in A 
\h* = on 

Recall that 5* is the normalized arclength measure on the circle Cn t of radius i?*. Up to now we have 



(2.75) 



proved, letting /3 in (2.72) 



/ - |V/i/| 2 < 47r 2 (AfiV) 2 (l + (l)) V / -\Vhl\ 2 +7rMNg 2 (R,)\loge\ + 

J A P i J A, P 



/MiVlog|loge| 
V e| log e| 



and arguing as in the proof of Lemma |2.3| we deduce 



1 

A P 



1 



V/i/l <Att z {MNY{1 + o{1)) I -|V/i*| +-KMNg*{R*)\\ogs\+0 

A P 



/MJVlog|log£:| 
V e|loge| 



(2.76) 



Compared to (2.431 there only remains to estimate the difference 

1 



V/i* 



A P 



\ |V^| 2 = / - \\7hJ 2 -h 
A 9 JaP 



where we have denoted h„ — hg t for short (see (1.291 for the definition of hg t ). 
Step 2. To estimate the above difference we note that 

-v(i(v^-v^))=v((i-^)v^) 

in A. Multiplying by h* — h*, integrating over A and using Cauchy-Schwarz, we obtain 

1/2 



- |V/i* - VhJ 2 < ( [ - \VK - V/iJ 2 
A P \JaP 



1 1 

A P \P 9 2 



2 \ 1/2 

\VhJ 2 \ . 



Using the definition of p and Lemma |2.1| it is straightforward to deduce that 



1 l - |2 

- \ VK - VftJ 

A P 



1/2 



< sup pg 



1 1 



.-i \ \P 9 
On the other hand, by Lemma |2.1| again 



1/2 



A 9 



1/2 



.4 9 



-|V/i„ 



1/2 



We conclude that 



^|VM 2 = / -|V^| 2 (l + o(l)). 
^ 9 2 JaP 



- |V/i*| 2 = I,(l + o(l)) 



and thus that (2.43) follows from (2.76) 



(2.77) 



□ 
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2.3 Interaction with the rotation field 



In this subsection we estimate the interaction of the vortices with the rotation potential (2.41 1: 

Proposition 2.3 (Interaction of the vortices with the rotation potential). 

Let v be defined in (2.11) and MN be as in (2.5). There holds 



2g 2 B ■ (iv, Vv) = 2nNF(R it ) + O 



e l/2|l ge|l/2 



O 



AH 



(2.78) 



where F is the potential function defined in (1.21). 



Note that F(R. lf ) is negative and proportional to e _1 so that the second and third term above are 
really remainders. 

As in the preceding subsection, the proof below will require many technical estimates. We thus find useful 
to give the core idea before going into the rigorous proof : Integrating by parts 



A 



2g z B ■ (iv, Vv) = / Fcurl(iv,Vv 



,H3 



F(iv, d T v). 



Now, F (see (1.27)) is defined in such a way that -F'(l) = 0(1). This has been proved in [CRYj to be a 
consequence of (1.22). Since on the other hand F(i?*) oc e , it is natural to neglect the boundary term 
above in a first approach. Then 



/ 2g 2 B ■ (iv, Vv) ~ / Fcurl(w, Vv) = V] / Fmr\(iv,V 

J A J A i ■ JB(p t j,2t) 



V) 



because |w| ~ 1 outside UijB(pij,2t), which implies cur\(iv,Vv) = there. Note that we neglect the 
thin boundary layer where \v\ is small for the sake of clarity. Then we use that t is very small to obtain 

/ Fcurl(w, Vv) ~ F(Pi,j) / cml(iv,Vv). 

Finally, using the fact that F is radial and by definition of the phase of v 



F(Vi,j) 



cw\(iv, Vv) ~ 2irF(R*) deg{v, dB( Ptd ,2t)} = 2irF(R* 



Collecting the above heuristics justifies (2.78), recalling that there are Af balls in the collection. 

That being said, let us go into details. We need to introduce a new potential function. Indeed, since 
the phase of our trial function is defined using the modified density p it is useful for the proof to define 
also a potential function accordingly. We denote 



F(r) := 2 
The function F satisfies 



1 



p(s) fls - ([fl] -u))-)ds = 2 / p(s)B(s) ■ e e ds 



R< 



and 



V ± F = 2pB, F(R < ) = 
-V Q VJ? ) = ~ 2V i B Cr) in A. 



(2.79) 
(2.80) 
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We can thus decompose F as follows : 

F = -Fin + F out 

where F n and F out are the solutions to the following elliptic problems 

f_V (±VF in ) = -2V (Be r ) in A 
{F in = on dA 

and 

-V(jVF out ) = 0iu.A 
' Fout = F(l) on OB 
F ont = on dB R . 

Both F{ n and F out are radial. The maximum principle implies 

|^out| < \F(l)\ 

and more precisely we have 

F out = F(i) r 



(2.81) 
(2.82) 

(2.83) 

(2.84) 
(2.85) 



where V is defined in (2.27). We sum up some useful properties of F in the following Lemma: 



Lemma 2.5 (Properties of F and F in ). 

Let F, Fj n and F out be defined as above and F in (1.27). The following properties hold 



VF in < 



< 


^|loge| 2 

L £ l/4 


(2.86) 


< 


L £ l/4 


(2.87) 


< 


£ l/4 


(2.88) 


< 


c 


(2.89) 


e 2 \ loge ' 



Proof. It is easy using the definitions of F and F to deduce (|2.86[ ) from (2.23). As F(l) = 
[CRYl Lemma_4.1]), (|2.87[ ) follows. A com binat ion of ( |2.86[ ), ( 2.84 ) and ( |2.87[ ) yields (|2J 
definitions of F and F out , Lemma 2.1 and (2.87) , we obtain 



0(1) (see 
Using the 



VF < 



C 



and 



VF out < 



e 2 | loge| 
C|loge| 



r5/4 



The estimate (2.891 then follows from (2.81 1 



□ 



We now present the 



Proof of Proposition 2.3 Step 1. We first replace g 2 by p in (2.41): Using \B\ < Ce 1 

C 



(g 2 -p)B- (iv, Vv) 



- I \g 2 -p\Z 2 \v< 

£ J A 



< 



C 



^ 2 |V0| : 



1/2 



v 2\ i/ 2 



(2.90) 
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From (2.42), (2.431 and (2.5) we have 



P em*<c(MNf+c^<c^L. 

A £ £ 



We also note that from the definition of p 



A P 



(9 2 -pY 



B.2 



where we have used Propositions 
bound to p respectively. We conclude 



^ln{r>_R} P 
and 



^Ce 1 / 4 / (g 2 -p) 2 <C£ 1/4 



B.4 



An{r>R} 

for the estimate of j An ^ r> ^ i (g 2 — p) 2 and the lower 



A 



g 2 B ■ (iv,Vv) = / pB ■ (iv.yv) + 



( (MN) 1 / 2 
\ en/8 



(2.91) 



We now compute, using (2.80), (2.81) and Stokes' formula, 



2 P B ■ (iv, Vu) = - / V ± (F in + F out )(iv,Vv) = / F in cuA(iv,Vv) - / V ± {F out )(iv,Vv). (2.92) 



(2.93) 



The second term is estimated as follows 



V ± (F ou t)(iv,Vv) = I V J -(f OMt )rV^= / F(1)^VT-Vh f . 



We then note that 



F(l) 



vr- v/i, 



A P 



But by definition (see (2.28)) 



< |^(1)| sup M 

A P i/2 



ivr 



(e i) 



1/2 



A P 



Vhtf 



1/2 



(2.94) 



sup^yy ^Csupp 1 / 2 < C{e\ logel)- 1 



A P 



-1/2 



A 



Also, £ < 1 with equality everywhere except in a boundary layer R < < r < R and in Ui,jB(pi.j , t). The 
former domain has an area of order s n with n large (we can choose it as large as needed) and the latter 
an area of order MNe 3 \ loge|. Thus 



/ (£ 2 — l) 2 < CMN£ 3 \ loge| 

J A 



Using also (2.87), Lemma 2.2 and Proposition 2.2 we obtain from (2.94) the estimate 



J A P 



vr • x/h 



< C(AfA0 3/2 £ 3 / 4 |log f 



(2.95) 



On the other hand 



F(i)-vr- vh f = F(i)[ / -vr-v/i/ 

P \JaP 



Ia \ IVlf J A P 



1 ivri 2 ] =o 



logej 

.1/4 
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because k = 



k = 0(1),F(1) = 0(1^) 



and 



-vr-v/i f = 

.4 P 



which follows from the definition of F and /i/ = on dA 



At this stage we have, gathering equations (2.91) to (2.95) and using (2.5) to estimate the remainders 
- J 2g 2 B ■ (iv, V«) = J F in cml(iv, Vu) + O ( ^^j^ ) • (2.96) 

SYep J2. To compute the remaining term we separate the contribution of the boundary layer from that 
of the bulk: 

/ F in cm\(iv, Vd) = / F in cur\(iv, V«) + / F in curl(w, Vu). (2.97) 

J A JAn{r<R} JAn{r>R} 

To estimate the boundary layer contribution we remark that for any r < R 

\F{r)\<^{r-R<)g\r) 
because \B\ < Ce^ 1 , p(r) = g 2 {r) for any r < R and g 2 is an increasing function. Also 

for any r < R, using the definition of T and the same arguments as above. We conclude that 

\F m (r)\ < -(r - i?<)5 2 M < Ce"" Y(r) 
e 



(2.98) 



(2.99) 



if r < i? (recall that i? = i?< +e" and that we are free to choose n as large as we want). We now compute: 



.An{r<_R} 



F in curl(w. Vu) 



<C I \F m \\Vv\ 2 

Ufl{r<fi} 



<C \F in \\W^\ 2 +C |F«HC| 2 |V0| 2 (2.100) 

JAn{r<R} JAn{r<R} 



For the first term we use (2.99) combined with the exponential smallncss of g in the boundary layer 



(Equation (B.4)). For the second one we use (2.99), (2.42) and (2.43) and conclude 



L 



AH{r<R} 



F in cwl(iv, Vv) 



< Ce r - 



(2.101) 



where it suffices to take n arbitrarily large in (2.16) to obtain that n' is arbitrarily large. 

For the bulk term we use that in A D {r > R} \ UijB(pij,2t), \v\ = c which is a constant, thus 

curl(iw, Vf ) = there. This yields 



/ Fi n cur\(iv, Vf ) = y / Fi n c\nl(iv, Vt> 

JAn{r>R} i:j JB( Piij ,2t) 



(2.102) 
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But, because of (2.18) and (2.891 



1,3 



y~] / {F in - F in (p itj )) curlO, Vv) 



< 



c 



e i/2| log£ |i/2^y fl(p . 3)2t) 



< Ce^llogel 1 / 2 / g 2 \Wv\ 
J A 

.MiVllogel 1 / 2 



- C- 



rl/2 



(2.103) 



Indeed, on B(pij, 2f) g 2 > C(e\ loge|) 1 and J A g 2 \^v\ 2 has been shown in the preceding subsections to 
be a O Finally 



Fin{Pi,j) / curlew, Vu) = 2irF in (p ii j) deg{v, dB(p id ,2t)} = 2TrF in (p itj ) 

JB(p itj ,2t) 

= 2irF(R*) + <D 



\oge[ 

£ l/4 



(2.104) 



by definition of the phase of v and ( 2.88 1 . The conclusion follows by gathering ( 2.96 ) , ( 2.97 1 and equations 
\2.1Q1\ to (12. 1041). □ 



3 Energy Lower Bound 



In this Section we provide the lo wer b ound announced in Theorem 1.3 The key to a lower bound matching 
the upper bound of Proposition 2.1 is the identification, in the energy of , 5 GP , of terms representing the 
kinetic energy of the vortices and their interaction with the rotation potential. This terms should then 
be bounded from below to show that our construction in the preceding Section is optimal. Ultimately 
we obtain the following 



Proposition 3.1 (Lower bound to the energy). 

Recall the decomposition 

\f GP = U ge im - u)B 



(3.1) 



valid in A and the definition of the reduced energy (2.35). Let I* be defined in 
function (1.26). Fore small enough there holds : 



1.31) and H be the cost 



> E<f u +£[u] - 0{e°°) > E% 



GP 



Ah 



-(1 + (1)). 



(3.2) 



The proof of this result will occupy the rest of the section. The main new ingredient with respect to 
[CRY] is a lower bound to the kinetic energy located outside the vortex balls, that is, relatively far from 
the vortex cores. In a first subsection we recall constructions from |CRY) and deduce some basic bounds 
that will be our starting point for the evaluation of the kinetic energy outside vortex balls. We refer 
to that paper for detailed explanations and comments. In a second subsection we show that the kinetic 
energy outside vortex balls can be bounded below using a variational problem related to electrostatics. 



We also present the analysis of the problem ( 1.31 1 as well as some important properties of the potential 
/i* associated to the minimizing measure 8* . We conclude the proof in a final subsection. 
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3.1 Preliminary Constructions 

The first inequality in ( |3.2[ ) has actually been proved in [CRY] Proposition 3.1] (note that this result 
stays available with our assumptions on O). The method is to combine the variational equation satisfied 
by g to obtain an energy decoupling, and the exponential smallness of ty GP in the complement of A. Our 
goal is now to bound £[u] from below to obtain the second inequality. 

To begin with, we define the energy 

T[v}:= ( 0*| V «| 2 + ^(l-|t;| a ) 2 (3.3) 

J A £ 

which plays a crucial role here. Indeed, a control on this energy allows to construct the vortex balls 
that are our main tool in this Section. The bounds that we use as starting point are the following (see 
Proposition 3.1 and Lemma 4.2 of |CRY| . again, this result stays valid under assumption (1.9 1, ( |1 . 10 1 
and 

Lemma 3.1 (First energy bounds). 

Let u be defined in (1.23). We hove 

T\u\ < Cs~ 2 (3.4) 
0{e°°)>£[u] > -Ce- 2 . (3.5) 

Our analysis uses the same covering of A by (almost rectangular) cells as in the upper bound Section, 
but we distinguish different types of cells : 

Definition 3.1 (Good and Bad Cells). 

We cover A with (almost rectangular) cells of side length C£|loge|, using a corresponding division of 
the angular variable as in Section 2.1. We note N oc -, the total number of cells and label the cells 

A n ,n € {1, ...,N}. Let < a < \ be a parameter to be chosen below. 
• We say that A n is an a- good cell if 

. 9 2 |V.| 2 + 4(1-H 2 ) 2 <^- Q (3.6) 



A„ 



We will denote the number of a- good cells and GS a the (good) set they cover. 
We say that An is an a-bad cell if 



x 5 2 |V.| 2 + f,(l-H 2 ) 2 >^e-« (3.7) 



We will denote the number of a-bad cells and BS a the (bad) set they cover. 

Note that the annulus A has a width i oc e| loge| (which implies that iV oc (e\ logej) -1 ) so that we are 
dividing it int o cel ls where there is much more energy than what would be expected from the localization 
of the bound (3.4) (namely ^ oc ^ lo ^ e ^ ) and regions of reasonably small energy. A first consequence of 
this is, using (3.4) and neglecting the good cells 



C C 
| log £r| " e 2 ~ e\ loge 

i.e. there are (relatively) very few a-bad cells 



N* < n^e^ < « N, (3.8) 
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The construction of vortex balls is feasible only in the regions of sufficient density. We thus introduce 
a reduced annulus 

A := {r : R> < r < 1} (3.9) 

with 

R> :=R h + ellogel -1 . (3.10) 
An important point is that from (B.6) and the definition (1.14) of p TF we have the lower bound 



9 2 {r) > 



C 



on A. 



e|loge| 3 

We can now recall the vortex balls construction (see CRY, Proposition 4.2] for the proof) 



(3.H) 



Proposition 3.2 (Vortex ball construction in the good set). 

Let < a < ~. There is a certain Eq so that, for e < £q there exists a finite collection {Bi}i e j :- 
{B(ai,ri)} i€l of disjoint balls with centers a% and radii fj such that 

1. [reGS a nA: \\u\-l\>\loge\- 1 } (lUUB,, 

2. for any a-good cell A n , ^2i b nA =£<z r i = ^ 

llogel- 5 . 

Setting di := deg{u,dBi\, if Bi C A fl GS a , and di — otherwise, we have the lower bounds 

,log|loge| 



J df g 2 \Vu\ 2 > 2tt Q - c?j \dk\g\ai) \\oge\ (l - 



(3.12) 



The second main tool that we need to introduce is the so-called Jacobian Estimate. For convenience 
we recall the result of [CRYl Proposition 4.3]: 

Proposition 3.3 (Jacobian estimate). 

Let < a < | and <fi be any piecewise-C 1 test function with compact support 



supp(</>) c AnGS a . 

Let {Bi} it -j := {B(ai, ri)} i€l be a collection of disjoint balls as in Proposition 
deg{u, dBi}, if Bi C A fl GS a , and di — otherwise, one has 



3.2 



Setting di 



^ 27rd i </)(a i 



i£l 



drift curl{iu, Vm) 



GS a nA 



< CIlVc 



e 2 \\oge\- 2 F[u]. 



(3.13) 



We refine our classification of cells 



Definition 3.2 (Pleasant and unpleasant cells). 

Recall the covering of the annulus A by cells A n , n £ {1, .., AT}. We say that A n is 

• an a-pleasant cell if A n and its two neighbors are good cells. We denote PS a the union of all 
a-pleasant cells and N F their number, 

• an a-unpleasant cell if either A n is a bad cell, or A n is a good cell but its two neighbors are bad 
cells. We denote UPS a the union of all a-unpleasant cells and their number, 

• an a-average cell if A n is a good cell but exactly one of its neighbors is not. We denote AS a the 
union of all a-average cells and iV A their number. 
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Remark that one obviously has, recalling (3. 



and 



N~ < 2N° < N. 



(3.14) 



(3.15) 

The average cells will play the role of transition layers between the pleasant set, where we will use the 
tools described above, and the unpleasant set, where we have little information and therefore have to rely 
on more basic estimates. 

We now introduce an azimuthal partition of unity that will allow us to avoid some boundary terms when 
integrating by parts in the sequel. Let us label UPS l a ,l € {1, . . . , L}, and PS 1 ™, m e {1, . . . , M}, the 
connected components of the a-unpleasant set and a-pleasant set respectively. We construct azimuthal 
positive functions, bounded independently of e, denoted by xY an( i Xm (the labels U and P stand for 
"pleasant set" and "unpleasant set") so that 



Xi 
X? 

, P 



E 



Xi = 



1 on UPS l al 

on PSl\ Vto G {1, . . 

1 on PS%, 

on UPS l a , VI € {1,.. 

1 on A. 



. ,M}, and on UPS l „, V/' + I, 
,L}, and on PS™' , Vm' ^ to, 



(3.16) 



It is important to note that each function so defined varies from to 1 in an average cell. A crucial 
consequence of this is that we can take functions satisfying 

C 



Vxl 



< 



s\ loge| 



< 



c 

£|log£ 



(3.17) 



because the side length of a cell is oc e\ log e| 
We will use the short-hand notation 



Xin ■ — 



Xout 



M 

/ j Am' 
rn — 1 



E 

Z=l 



(3.18) 
(3.19) 



The subscripts 'in' and 'out' refer to 'in the pleasant set' and 'out of the pleasant set' respectively. 

We want to use the Jacobian estimate of Proposition |3 ,3| with <fi = XinF, which does not vanish on dB 
and has its support included in A which is larger than A. We will need one more construction to make 
this possible : We introduce two radii R* ut and R~ ut as 



R c 
Rf 



l-ellogel" 1 , 
R> +e\ logel" 1 . 



(3.20) 
(3.21) 



Let £in( r ') and ^ ut('') be two positive radial functions satisfying 

Cout(f) 

Sout(r) 

6n + 6out = 1 on A 



1 for R- Ut < r < R+ Ut , 

for i?< < r < i?> and for r = 1, 

1 for i?< < r < R>, 
for R mt < r < i?+ t , 



(3.22) 
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Moreover, because of (3.20) and (3.211, we can impose 



|V6n| < 



C|loge| 



|VUt| < 



C|loge| 



e e 
The subscripts 'in' and 'out' refer to 'inside A' and 'outside of A' respectively. 



(3.23) 



In the sequel {Bi} ieI := {B(ai,ri)} i£l is a collection of disjoint balls as in Proposition 3.2 For the 
sake of simplicity we label Bj, j E J c /, the balls such that Bj C A D GS a . 



An important step towards a lower bound to £[u] is to remark that, integrating by parts, 
- / 2g 2 B ■ (m,Vu) = / Fcurl(ra, Vu) - F(l) / {iu,d T u) 



A 



A 



(3.24) 



and to note that the boundary term above can be neglected. In |CRY) (see equation (4.100) therein) we 
have proved 

<c(\lo g en£[u}^ + ^l^nu]^ 



F(l)(iu, d T u) 
and thus we deduce from Lemma 13 . 1 1 that 



F(l)(iu,d T u) 



< C 



loge|V4 



r5/4 



(3.25) 



Note that this term is much smaller (in absolute value) than the lower bound we are aiming at. 
Our first lower bound is the intermediate result (4.86) in [CRY] . that we reorganize to obtain : 



/ g 2 \Vu\ 2 + FcuY\(iu.yu) > (I-27) / 6n£ 2 |V U | 2 +J2 ^ 2 (% 



|d J |(l-27)(i-a)^(a i )|log £ |(l-C 



,l0gjl0g£| 

|loge| 



lXout(o,-)|dj||loge|(l-o(l)) 



iljCuA" ,)F(<ij) 



I &?Ut9 2 \^u\ 2 ~ ( drUt\F\\Vu\ 2 -C f da\F(l)\\(iu,d T u)\ 

J A JA JdB 

Arg 2 \u\ 2 - C\ lDge| -1 .F[u]. (3.26) 



+ h-6) J^dfg 2 \Vu\ 2 -^ 



UPS a UAS a 



We emphasize that we have kept the kinetic energy contained outside the vortex balls (first term in the 
right-hand side, on the first line) that was neglected in [CRY]- The parameters in ( |3.26[ ) are chosen as 
follows : 

log | log e | .log | log e| 

a = a- 



7 = 25 



(3.27) 



|loge| | log e | 

where a is a large enough constant (see belo w). 

Using the estimate on the number of bad cells (3.8| and the upper bound on g 2 \ u | 2 = |\I> GP | 2 of Proposition 
IB. II we have 

1 

Also, using [CRYl Equation (4.22)], 



? 2 M 2 <c- |loge|£Q 



UPS a UAS a 



e 2 log | log £ 



7 /yiv„|*-/^.lW> C ™//|v«f 
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Plugging these estimates and (3.25) in (3.261, using Lemma 3.1 to estimate the last term, and taking 
a > 2 we obtain 



/ g 2 |V M | 2 +Fcurl(m, Vu) > ( 1 - C l °f) l ° g f \ f ^ n g 2 \Vu\ 2 +C Y i si > (a j )x aa t(.a j )\dj\\ loge 



^27rxin(a J ) 



AXUjejBj 

M,l^ 2 K)|log £ | (l-C ^y 1 ) + d j ^(a j )F(a j ) 



-C 



loge| 



e 2 log | log e | 



(3.28) 



The critical speed 2(3ne 2 \ loge|) 1 is (roughly speaking) defined as the first speed at which the terms 
on the second line of (3.28) all become positive. This corresponds to the speed at which the vortices 



cease to be energetically favorable. In |CRY| we were above the critical speed, so the preceding lower 
bound was enough for our purpose, because the terms on the second line were positive. Here we are in 
the opposite situation where the vortices can become energetically favorable if they are suitably located 
in the annulus. Thus, a lower bound to £[u] requires an upper bound to the number of vortices (more 
precisely, to the sum of their degrees). We will provide this upper bound in the sequel. 

As is standard in such problems, we need to distinguish between different types of vortex balls. First 
we need to distinguish the vortices lying close to the inner boundary of A from those in the bulk. The 
vortices close to the inner boundary of A have positive energy and can thus be neglected in the lower 
bound. However their energetic cost is not large enough to show that there are few vortices of this type. 
Also, vortices in the bulk can be energetically favorable only if their degrees are positive. On the other 
hand, in the range of S7 that we consider even a vortex of positive degree can lower the energy only if it 
is close to the circle Cr, . 
Thus we first divide J into two subsets 



J in = | j G J , K| > i?buik = Rh + e\ loge|fii /2 1 

-^out J \ Jin- 



We next divide J; n into three subsets 



J- = {j € J in , dj < 0} 

J+ = jj G Jin, dj > and ||o 3 -| — R*\ > e|loge|^y 4 | 
J* = J in \J-\J+. 



(3.29) 
(3.30) 

(3.31) 
(3.32) 
(3.33) 



With these definition we can state our lower bound to the kinetic energy contained 'far' from the vortex 
cores : 

Proposition 3.4 (Lower bound to the kinetic energy outside vortex balls). 

Recall the definition of I* (1.31). Let J* be as above. There holds 



^g 2 \Vu\ 2 > L^XinKK-) f >( l + 0(<))- C g^g^ 



(3.34) 



The energy evaluated in Proposition 3.4 is associated with the supcrfluid currents that the vortices 
induce in the condensate. As is well known the current generated by vortices of negative degrees 
could compensate that generated by vortices of positive degrees. In the sequel we prove that there are 



relatively few vortices of negative degrees (see (3.42) below). We will thus deduce that this phenomenon 



does not affect the leading order of the energy in the situation we consider. 
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Note also that the minimization problem appearing in our lower bound involves vorticity measures with 
support on Cr, . Such a restriction of the set of admissible measures will be proved to be favorable 
because most vortices have to be located close to the circle Cr, where they are energetically favorable 



(see (3.43) below) 



We begin the proof of Proposition |3.4| with the following lemma that gives lower bounds to the 
energetic cost of the vortices, depending on their degrees and locations in the annulus. The proof is 
postponed to Appendix A. 



Lemma 3.2 (Energetic cost of the different types of vortices). 

For any j £ J out 



1 



\d 3 \-g'(a 3 ) \loge\ 1-C 



Y log|loge| 
|loge| 



For any j £ J h - 



1 



2^3 \y v*j 

Moreover, if j £ J_ U J + 



Recalling that 



djtfia-) |loge| + d j ^(a j )F(a j ) > \d 3 \H(R*)(l + > -C\d 3 \ 



Qi 



1/2 



-\d 3 \g z (a 3 ) |loge| + d 3 ^ n (a 3 )F(a 3 ) > \d 3 \- 



(3.35) 



(3.36) 



(3.37) 



log I log g I 
|loge| 



< Oi < 1 



we can use Lemma 3.2 to simplify further (3.281 



f g 2 |V M | 2 +FcurlKV M ) > (l- l0 Sl l0 g £ h f ^ n g 2 \Wu\ 2 + Cj29 2 (a 3 ) X out(a 3 )\d 3 \\loge\ 

.1/2 



Cfl ' 

2*H{R.) (1 + o(l)) £ Xin(a 3 )d 3 + Y, Xin(a 3 )\d 3 



C 



jeJ, 



je.J+u.J- 



loge|e a 
e 2 log I log el 



(3.38) 



Note that by definition of £ in , £ in (a 3 ) = 1 for any j £ J*. Adding f A |y (l — |u| 2 ) — F(l) J dB (iu, d T u) to 
both sides of (3.38) and using (3.25) we obtain 



E[u] > 1 



log I logs I 
|loge| 



(aj)Xout(aj)\d 3 \\loge\ 



2nH(R*)(l + o(l))Yx i n(a 3 )d 3 + CD ^ £ Xin^OMjl ~ C J^li ■ ( 3 - 39 ) 
jeJ* jeJ+uj- -Logiiogti 



On the other hand, combining the upper bound to the GP energy of Proposition |2.1| and the first inequality 

(3.40) 



in (3.2 ) we have 



n 2 

S[u] < -C% 



We deduce from the above 



^ Sli \£ £ 2 log | log e | y e 



(3.41) 
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The second inequality holds true if we choose a > 3, which we now do. From (3.39) and (3.40) we also 
deduce (note that the last term in ( |3.39[ ) is a remainder because of (3.41 1) 



jeJ- jeJ, 

Y Xin(aj)\dj\ < Ctt\ /2 ^2 Xin(Oj)|dj 

jeJ+ jeJ, 



(3.42) 
(3.43) 



One can interpret (3.41), (3.42) and (3.43) as follows (recall that fix < 1): There are at least C^ 1 
essential vortices in A (meaning vortices with nonzero degree). Most of them are of positive degree and 
close to the circle Cu t . 



3.2 The electrostatic problem 



A difficulty in the proof below is to define a vorticity with support on Cr^ starting from u. A possible 
track, following [ABM], would be to use the Jacobian estimate and a first rough upper bound to the 
number of vortices to obtain some compactness for the vorticity measure of u. Using (3.43) one would 



then show that the limit measure has its support on the circle of interest and obtain the lower bound 
by a lower semi-continuity argument. Such a strategy is difficult to adapt to our setting because the 
geometry of our domain strongly depends on e. To obtain non trivial limits one should rescale the 
annulus to work on fixed domains. Also the weights appearing in the energy would complicate the 
argument (recall that g 2 is very small close to dBn < ) . 

We follow another route that does not require any compactness argument. In this subsection we prove 



that the left-hand side of (3.34) can be bounded below using the electrostatic energy of a suitable 



modification of the vorticity of u. The method that we use to conclude the proof of Proposition |3.4 
in Subsection 3.3 requires precise informations on the potential h, = hg t generated by 8, according to 



(1.29). We provide this information in Proposition 3.5 below 



Let us describe the electrostatic energy that will serve as intermediate lower bound. For technical 
reasons it is necessary to reduce the domain on which we work to *4b u ik where the density is large enough. 
For any Radon measure v supported in .4buik we define h v as the unique solution to the elliptic problem 



-V iV/iJ = v in Abu 



lk 



h v = on cl4bu 



(3.44) 



lk- 



and introduce 

I{v) := [ \ \VhJ 2 . (3.45) 
We will later be interested in the minimization problem 

I, := inf I(v) = 1(8,) (3.46) 

i/€P,J i/=l 



where the infimum is taken over the set T>, of positive Radon measures with support on the circle of 
radius i?*. The fact that 8, is the unique solution of (3.46) will be proven below (Proposition 3.5). 
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The definition of the measure whose energy will be used in the lower bound goes as follows : Let j 
be the modified superfluid current 

~. _ j (iu, Vu) in _4 bu ik H GS a \ Uj eJ Bj ^ 
[ otherwise . 

The measure fi is the vorticity associated to j 

A := curl(j). (3.48) 
The following lemma is a key ingredient in our analysis : 

Lemma 3.3 (Lower bound via an electrostatic energy). 

With the above definitions, there holds for e small enough : 



I ^\Vu\ 2 >{l-C\\o E e\- x )I{Jl). 



Proof of Lemma \3.3[ We need to define T as the solution to 



(3.49) 



-v (^vrj = o in Auik 

r = on dB Rbalk < :! - r > l) > 
v T = 1 on dB. 

Explicitly : 

g 2 (s)s~ 1 ds 

r( f) = J ^ b - k . (3.5i) 

jfl. bulk 5 2 ( s ) s lds 

Note that this function is not the same as that defined in (2.27), which is not confusing because the latter 
will not be used again in the sequel. We recall that in ^4buik H GS a \ Uj e jBj, £i n = 1 and \u\ is close to 
1 according to Item 1 in Proposition |3.2| Thus 



6nff 2 |Vu| 2 > / g 2 \Vu\ 2 



/ 

> (l-qioger 1 ) / g 2 H 2 |v u | 2 

> (1-qioger 1 ) / g 2 \(iu,Vu)\ 2 

JAb^nGS^UjejBj 

= (1-qioger 1 ) f g 2 \]\ 2 . (3.52) 



We set 



By definition 



curl (J+^V X /,, ) = in // 'Uu.nJ 

i+iv^/^j-r = 0. (3.54) 

98 V 9 
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Hence there exists / G -ff 1 (^4buik) such that 

3 = 



Using the fact that fa is constant on the boundary of *4buik, we have J A ^ V^/^ ■ V/ = and thus 



-Abull 



Next we note that 



9 2 \3\ 2 > 



1 



1 



•Abuik 9 



2 |V/ A | 2 . 



2 v/i p -vr 



Mbuik 5 

because hp, = on cMbuik = dBR buik U 9B. We thus have 

1 



-Abuii 



-Abuik 9 



r|V/ A l 



> 



1 



-Abuik 5 



Combining with (3.52 1 we conclude that (3.49) holds. 



□ 



Our next task in this subsection is to give some details on the minimization problems (1.311 and 



(3.46). In particular we prove that the minimizing measure is in both cases the normalized arclength 
measure on Cr^ and we compute the associated potential explicitly. We also show that considering the 
problem on the reduced annulus *4buik does not change significantly the energy. 

The following proposition contains probably only facts known from potential theory. Indeed, apart from 
the weight g~ 2 , the minimization problems we are considering fall in the general context of |ST] (see e.g. 
Theorem II. 5. 12 therein and the discussion in [ABM1 Example 5.1]). We nevertheless provide a short 
proof for the sake of completeness. 

Proposition 3.5 (The electrostatic problems). 



1. The minimization problems (1.31) and (3.46) both admit the normalized arclength measure on 
Cr, (denoted 8*) for unique solution. Moreover there exists two positive constants C\ and Ci 
independent of e such that 

Ci<h<C 2 . (3.55) 



2. We denote h* — hg t and h* — h$ t the fields associated to <5* by equations (1.29) and (3.4-4). Both 
are radial piecewise C 1 functions. Explicitly we have 



(f R ^{s)s-Hs 



Jr<9 2 (s)s l d& 
Jr9 2 (s)s~ 1 ds 



h*(r) = h*(R*) < 



for R < < r < R* 



(3.56) 



ds 



and 



for R* <r < 1 



h*{r) — h*{R*) < 



-Rbuik 
^,2 



'-(s)s- 1 ds 



for i?buik <r < R* 



J r 9 2 (s)s x ds 
Jr, 9 2 (s)s~ 1 ds 



(3.57) 



for R* < r < 1 
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where 

ft* (-R*) 
ft* (-R*) 

5. Tfte following estimate holds 



2^< 3 2 (^' 



L d.s 



2tX 



,, q 2 (s)s 1 ds 



1/2 



(3.58) 



(3.59) 



(3.60) 



Proof. The results for the problem (3.461 are exactly similar to those for (1.31), we thus only prove the 
later. 

Let us start with the general problem of minimizing I(y) over the set T>e of Radon measures supported 
on £, a compact subset of A : 

I E = inf I{v). (3.61) 

veV E ,f u=l 

The existence of a minimizer ve to such a problem is classical (see e.g. |ST| ). We denote He the function 
associated to v E by formula ( 1.29 ). Next, computing the first variation of /, we observe that there holds, 
for any v G T>e 

' ' (3.62) 

Je Je 

ier. We deduce that He is constant on E. Using the above equati< 

(1.29) with v = ve we see that 



/ vh E = 
Je 

where As is a Lagrange multiplier. We deduce that He is constant on E. Using the above equation and 

1 



l\E 



A, 



Vft, 



(3.63) 



We turn to the proof of Item 1. Let ft* be a solution to the minimization problem. Using (3.63) and the 
definition of ft* we have 

-V(p-Vft*) = in AnB Rst 
ft, = on 9B_r < 
ft* = /* on dB Rt . 

This implies that ft* is radial on A H B Rt . A similar argument yields that ft* must be radial also on 
A \ Bfj r and thus radial in A. Then the associated minimizing measure vc R is radial also. We conclude 
that <5* is the unique solution to the problem (1.31). The bounds (3.55) are proved by noting that, using 
C29l 



^|Vft*| 2 
A 9 



sup 

<pecl(A) 



Ia ■ v ^> 



1 



, „ 9 sup 
4w 2 Ri 4, e ci(A) J 



IA g 2 I v VI 1,1 "* pfcOJlAJ J^t ff 2 

The upper bound follows because g 2 < C(e| log e|) — 1 on A whose thickness is of order e\ loge|. A small 
computation shows that the above supremum is bounded. The lower bound is proved via a trial function 
for the maximization problem. For example one can use the test function used in Remark |1.1| Item 2. 
Proving that the function defined in the right-hand side of (3.56) and (3.58) is a solution to ( |1 .29 1 with 
v = 5* is a straightforward computation. By uniqueness we deduce that it must be equal to ft* . 
Finally (3.60) follows from (3.55) and the formulas (3.58) and (3.59) once one has noted that 

l-Rbulk — R<\ 



g 2 (s)s- L ds < C 



e|loge| 



1/2 



□ 
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3.3 Completion of the proofs of Propositions 3.1 and 3.4 

In what follows we denote 

D := 2tt ^2 Xin{dj)dj 



and 

We also decompose ji as 



^ := D5*. 
/i = /i* + fi = DS* + (l. 



It is useful to recall that from (3.41) (3.42) and (3.43) we have 



X 

e 



(3.64) 

(3.65) 
(3.66) 

(3.67) 



We are now ready to finish the 



Proof of Proposition \3.4\ Starting from Lemma 3.3 there remains to evaluate I(fi). The two main ingre 
dicnts will be the Jacobian Estimate and the properties of h* proved in Proposition 
We have 



3.5 



1 



alburn 9 



1 



Now, 



^ibuik 9 
1 



1 



-Abuik 9 



2 Vhn- V/i M . 



1 



Abuik 9 



-Abuik 9 



VV| 2 =D 2 l(5*) = D' 2 l4l + <D(tl{ /2 )) 



(3.68) 



(3.69) 



by definition and use of (3.60). To obtain an appropriate lower bound it is thus sufficient to compute the 

ve 

Kfi-DK{R*) 



second term in the right-hand side of (3.68). Using (3.44) and Proposition 3.5 we have 
1 



-Abuik 9 



h^fi = D 



-Abuii 



h*ji — D 



-Abuii 



-Abuii 



/i*<5* = D 



-Abun 



(3.70) 



We now compute from the definition (3.47) of j 



-Abuu 



-4 b uii 



GS a 



(iu, Vu) ■ V 1 - (xinK) - / (in, Vu) ■ V 1 - (xoutK) 



E 



(iu, Vu) ■ V ± h* 



/ VXinh - [ (iu,Vu) ■ V x (xouth) +J2 I K Vu) ■ V^k. (3.71) 

JGS a JgSc, - c T^B.-nylbuik 



Let us first show how to estimate the second term. The integral is actually located on 

-Abu* n (AS a u (e/ps q n GS* Q )) 
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(see the definitions of x out and h*). Recalling that g 2 > Ce 1 | log e| 3 and |Vx out | < C (e| log e|) 1 there 
we have 



GS a 



(iu,Vu) ■ V 1 - (xoutK) 



, c[ IMki + || V ^|| \ f | u || Vu | 

£ l l0 S £ l / JA hv i k n(AS a u(UPS a nGS a )) 



< C\\oge\ 2 [ 

J A 

< C\\oge\ 2 V [ 

J A 



AS a U(UPS a nGS a 



g 2 \u\\Vu\ 



9 2 \u\ 2 



AS a U(UPS a nGS a ) 



+ 



C|l0g£| ; 



V JAS a u(UPS a nGS a 



g 2 \Vu[ 



(3.72) 



where 77 is a parameter that we fix below and we have used the explicit formula (3.571 for the bounds on 
h„ and \\7hJ. We recall the basic estimate 



9 2 \Vu\ 2 < 



AS a U(UPS a nGS a ) 



C 
72 



coming from (3.4). Also, using g 2 |«| 2 < C(e\ loge|) 1 



AS a U(UPS a nGS a ) 



g 2 \u\ 2 <Ce a , 



as a consequence of (3.15) and (3.14). Choosing 



rj = e 



-l-a/2 



and using (3.27) we obtain 



GS a 



(iu,Vu) ■ V- 1 (xout^*) 



< C 



I loge| 2 

£ |l ge|5/2" 



(3.73) 



For the third term in (3.71) we use exactly the same kind of argument. We now use the smallness of 



the set covered by the vortex balls instead of the smallness of AS a . Indeed, using the division of A into 



N oc (e| log e| ) 1 cells and Item 2 of Proposition 3.2 

N 



N 



n=l j.BjnA n ^0 



(3.74) 



n=l 

_2|„,|2 ^ n(A l„„-l\-l^ 



It follows that (recall that g 2 \u\ 2 < C{e\ loge|)~ 

g 2 \u\ 2 <C\\oge\~ 12 . 



Using this fact and arguing as in (3.72) above we obtain 



I (iu,Wu) ■ W^h, 



jeJ" B 3 



< C 



e\ loge\ 



(3.75) 
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The first term in (3.711 is esti mated using the Jacobian Estimate : The function Xinh* satisfies the 
assumptions of Proposi tion 3.3 In particular, its support is i nclud ed in _4b u ik fl GS a C A f~l GS a . We 



thus have, using again (3.57) to bound the gradient of h„ and (3.4) 



/ MXin^* = 2tt V" Xin(ajKM a i) + ( -jj— 
JGS a £[J in V e l io g 



(3.76) 



Provided a is large enough (which we are free to decide) we thus finally have, inserting (3.73), (3.75) and 
(|3776b in pTfTl) 



I hp, = 2n^2 Xin{a.j)d K{a 3 ) + O ( J — 

■Mbulk j'GJin \ fc l iU 6 fc l 



(3.77) 



Then, going back to (3.70) 



J-A buik 9 \ je.hn I \ E P § 

On the other hand by definition of J* and the explicit formula for h* 

\h*(dj) — h*(R*)\ < CQ 1 ^ 4 for any j e J*. 
Recalling that h* (i?* ) = max_4 bulk h* we have for any j e J_ 

rfj (h*(a,j) — > 0. 

Also, for any j 6 J + , 



(3.78) 



It follows that 



> -CD £>ftj 



)V4 



el log el 3 



+ 0( £ - 1 |log £ |- 3 ) 



(3.79) 



where we have used (3.67). Going back to (3.68) and (3.69), combining with (3.49) and (3.67), the result 
is proved. 

□ 



With the result of Proposition 3.4 in hand it is an easy task to complete the 
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Proof of Proposition\3.1\ Collecting (3.34 1 and (3.39) we have 



)l/4 



s[u] > ^2tt |^ XinMdjj L(i-cn{ 

+ 2icH(R*)(l + o(l))J2xU*j)d j 

jeJ, 

( X, ./. V' Ul r d J c |loge|e" 
£ | log e | 3 e 2 log|loge| 



+ 2wH(R*) (1 + o(l)) XinKK" - C 



log | log £ | 
£ 2| log£ |l+5- 



(3.80) 



We have used ( |1.11[ ) and (3.27) for the second inequality. There only remains to minimize the sum of the 
first two terms above with respect to 27r^V gJ Xin{a>j)dj and use (1.10) and (1.11) to obtain 



(3.81) 



The second inequality holds true because we choose a large enough and H(R^) oc — Vt\E~ x , This concludes 

the proof of the second inequality in (3.2). 

Note for later use that in (3.80) we have dropped the term 



Keeping this term and combining the lower bound we have just proved with the upper bound to the 
energy we obtain the estimate 



cY J 9 1 {a j ) X o^{a 2 )\d 2 \\\oge\ < 



(3.82) 



□ 



4 Asymptotics for the vorticity 



In this Section we provide the proofs of Theorems |1.1| and |1.2| Actually, most of the ingredients are 
scattered in the preceding subsections. In particular it is rather straightforward from the proofs in 
Subsection 3.3 to deduce 

Proposition 4.1 (Asymptotics for the modified vorticity). 

Let jl be defined as in (3.4-8). We have for any test function 4> € C^(^lbuik) 



H(R* 
2L 



/ <5*0 










Mbuik 







1/2 



(4.1) 



uniformly in <f>. 
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Proof. We first remark that, combining (JlTTJ), (|2T2j) , (|372j> and ( |3.80 ) 

2 



41* 



-(1 + 0(1)) > 



2tt ^ Xin(aj)^ J /* (1 



o(l)) 



•0(1)) ^ Xin(Oj)dj -O 



of 



We deduce that 



2tt 



Xin(a,K = -ffl^(l + o(l)) 



and thus, going back to (3.67) 



D = -^(1 + 0(1)). 



(4.2) 



(4.3) 



(4.4) 



Next we note that, up to now, we have neglected one useful term in the proof of the lower bound, namely 
the third term in (3.68). Keeping this term and using (4.3), Equation (3.80) reads 

H(R, "> 2 



AI, 



"(1 + 0(1)) 



1 



9* 



\Vhp\'(l + o(l))-o 



£[u] > 

Thus (recall that fx = ft — D6*), using (2.2) and the first inequality in (3.2) 



(4.5) 



1 



-Abuik 9 



O 2 



which implies the result via (3.44) and (4.4). 



□ 



Recall that fl is defined by neglecting the superfluid current in a region where no convenient lower 
bound to |u| is available and considering the vorticity associated to the remaining current. One can see 
this procedure as a regularization of //, because we actually neglect the regions where we expect the phase 
of u to be singular (therefore |u| to go to zero). Proposition 4.1 states that, modulo this regularization, 
one can estimate the vorticity in the dual norm of 



1/2 



This is the natural norm associated with the minimization problem (1.31) defining the contribution of 
the vortices to the energy. Note that such a control is not exactly stronger than a control in (C^^buik))* 
norm because of the strong inhomogeneity of the weight g~ 2 (see Remark 1.1). It would be stronger for 
example if g 2 was a constant because of the embedding (in two dimensions) of iJ _1 (^lbuik) in (C<~ C^bulk))* • 
The necessity to use the L°° norm of the gradient of test functions to bound the left-hand side of ( |1 .34 ) 
appears when approximating \i by fl in the proof below (that is when justifying (1.54) rigorously). The 
main task is now to estimate the effect of the currents that we have neglected when defining fl. 

Proof of Theorem \l.l\ Let <p £ C^(Ahuik)- We compute 



fJ0 



- / (iu,Vu)V (Xin^ + Xout<; 



■I 



(iu, V^V 1 (xi„0) 



A 



-Abulk\U j^jBj 

(m, (xout<i 



(iU,Vu)V X (Xin^) 



-AbuikHUje,/ Bj 



(4.6) 
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The first term is estimated using Proposition 4.1 (note that Xin<f> nas its support included in -4buik)- By 
definition 



A b u 1 k \ -B S a \ U j £ J B j 

_ H(R,) 
2L 



(m, Vu)V ± (Xin0) = / MXin^ = - 



27, 



<5*Xin0+O 









''-Abuik 





1 



•Abuik 9 



1/2 



27* 



+ 



-4 bu lk f 



■|V# 



-4-buii 

• | A.( ! 1, 

1/2 



1 

^buik y 



1/2 



Ii°°(^bulk)" 



•0(0x1 log ff|||V^|U ro( ^ lk) ). (4.7) 



We have used the fact that Xin = 1 on Cr„ n PS a . To pass to the second line it is thus sufficient to note 
that the length of Cr, \ PS a = Cr, H AS a is of order N® e\ loge|. Also, we have 

HWl^(A) < CellogelUV^IIioc^) 



because |.4| oc e\ loge| and <j> vanishes on OA. Recalling (3.8) we obtain the third line of (4.7 1. 

We now show that the other terms in (4.6) are remainders, arguing as when dealing with (3.71) in the 

proof of Proposition |3.4| 



The second term in the right-hand side is estimated exactly as the sum of the second and the third term 
in the right-hand side of (3.71 1, using the small area of the region covered by the vortex balls. The result 
is 



(w,Vu)V x (x 



<C|loger 3 ||V0|| io 



■(^bulk)- 



Finally, the last term in (4.6) is estimated exactly as in (|3.72|), (|3.73|). We obtain 

/ (i^VttJV 1 - (Xou 

Ja 



logs 



i/2' 



which concludes the proof, taking a large enough. 



(4.8) 



□ 



We conclude this section by the 
Proof of Theorem The vortex balls entering in the definition of /i e are those defined in Proposition 



3.2 from which we discard the balls that are not included in >4buik, that is the balls labeled by j G Ji T: 



(see Equation ( |3.29| )): 

He = 2ird.j5 aj . 

For the statement of the Theorem we have renamed J- m = K. 

Let us first note that, using Theor em [H] ( fL"40| ) is a consequence of (II. 391), we thus only prove the latter. 
From the Jacobian estimate (3.13) and Lemma 3.1 we have, for any <fi £ C^(A.hu\k H GS a ) 



N 2-Kolj(j)(ai) — I d 

jGJ iD JGS a nA 



r(f> curl(iu, Vit) 



< CHVc 



I log e | 



where C does not depend on <j>. To conclude the proof we only have to extend such a statement to all 
test functions (f> € C\ (^4buik H GS a ) . Let us pick such a function and write 

(jl - fl e ) (j) = / (fj, - H e ) Xm<t> + / MXout^- d jXout(aj)4>( a j)- 

bulk ^-^bulk ^-4-bulk jGzJin 
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For the first term we can use the Jacobian estimate because Xin<P has support in GS a . This yields (using 
also Lemma 3.1 1 



(fl - /i e )Xin</> 



The second term is 



< C|loger 2 ||V0|| Lo 
(iu, Vu)V J -(xout<> 



and has already been estimated, see (4.8 1. For the third term we recall that on -4tmik 

est 1 ' 2 



9 2 ~ P TF > 



e\\oge\ 



thus, using \<f>(aj)\ < Ce|loge| ||V0|| Loo and (3.82 1 



d jXout(aj)<f>(aj 



< Ce^logel^^/^^KOxoutKOI^-l « n? /2 |loge| ||V0|| 



This concludes the proof. 



□ 



Appendix A : The Cost Function and the Vortex Energy 



In this appendix we study the cost function (1.261 and provide the proof of Lemma 3.2 
It is convenient to define 



,TF 



and study the related function 



H TF (r) := || log e\p TF (r) + F TF (r) 



where 



F TF (r) := 2 f ds B^{r) ■ e> TF (r) = e 2 Sl 2 f ds [{Is - ([SI] - lo tf ) s" 1 ] (s 2 - Bfc). 



(A.l) 
(A.2) 
(A.3) 



In order to investigate the behavior of the infimum of _ff TF inside the bulk, it is convenient to rescale the 
quantities and set 

z := eil(r 2 - Rl), (A.4) 

so that z varies on a scale of order one, i.e., more precisely z € [0, 2/t/tt\. With such a choice the gain 
function can be easily estimated: 



F TF (r) = 



2 si 2 r 2 - R i 



dt t [fi(t + Rl) - [SI] + uj tf ] (t + Rl 



2 n 2 r 2 - R * 



dt t 1 Sit 



+ 0(1) 1 



f-KeSl 



+ t 



1 

2eJ 



+ 0{e) 1- 



^JlteSt eSl 



1 

Ye 



— ) 



•0(|loge|) = 



■O(|loge|), (A.5) 
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where we have used the approximation [1 — 0((efl) x )] 1 = 1 + 0((efl) 1 ). 
Applying the same rescaling to the energy cost function, we thus obtain 



12e ' 



where (recall that ft = (2 (3tt) 1 - Ct^j e~ 2 \ logel" 1 ) 



H TF (z) = z I - -30i -2z 

7T 



0(e|loge|). 



(A.6) 



(A.7) 



Let us denote, for z e 



fc(z) := 



3fii - 2z 



(A.8) 



It is straightforward to see that this function takes two local maxima and two local minima in 
The maxima are at -?= and at 

2 



zi = 



3v^ 



The minima are at z = with fc(0) = and at 



zi = 



3v^ 



1 fii 

V — . 

9tt 2 



1 fti 
9^ + ~2~' 



Computing this local minimum we obtain that k(z 2 ) < if and only if fix > 0, thus k(z 2 ) is the absolute 
minimum in this regime. More precisely, for |f2i| <C 1 we have 



k(z 2 ) = — + o(nj) 



and thus, defining i?* > by 
we obtain, for (e\loge\) 1/2 < jOi j < 1 

tf TF (i?*) 

Also 



Rl := R 2 h + (sty- 1 z 2 



and 



fc"(* 2 ) 



fc'(0) 



1 



(A.9) 
(A.10) 

(A.ll) 



•O(fii) 



3fii. 



Recalling that fc(z) increases from to zi and from z 2 to 1 and decreases from z\ to z 2 , we have for any 
z and an appropriate choice of constants fci, fc 2 , £3 such that |z — z 2 | > fciSl^ 4 and \z\ > /c 2 f2^ 2 

fc(z) > k z n\ /2 . 



Collecting the preceding facts we have proved 
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Proposition A.l (TF vortex energy). 

Let fl be of the form 

2 

37re 2 |loge| e 2 |loge| 
with (e| log e|) 1/2 < O x < 1. Let f e A TF . We have 



H rt (r) > 



4 v / 7re 

n 1/2 



O 



H Li! (r) > C~-^ifr>R h + Cs\\oge\n\ /2 and\r-RA>Ce\\ogs\n[ 

£ 



1/4 



(A.12) 



(A.13) 



(A.14) 



We now compare the cost TF function H TF to the original function appearing in our analysis, H. 
The following result is a part of the proof of Proposition A. 2 in |CRY) and is sufficient for our purpose. 

Proposition A. 2 (Comparison of the cost functions). 

Let ijj and g = qa.ui be defined as in Proposition \l . 1\ Let ft be as above and H be the cost function defined 
in (1.26). For any r € A (see definition {3.9))) there holds 



\ 9 2 (r)-p TF (r)\ < 
\F(r) - F TF (r)\ < 



C|loge| 5 / 2 



£ l/2 
C 



e\\oge\ 



and thus 



\H{r)- H lt (r)\ < 



C 



e|loge| 



(A.15) 

(A.16) 
(A.17) 

(A.18) 



We are now equipped to present the 



Proof of Lemma 3.2 The proof of (3.35) is based on a simple estimate from [CRY} Lemma 4.1] that we 
recall 



\F(r)\ < Cmin 



\r-R 



*r(r),i + 



C\r-l\ 
e 2 1 log el 



Note that this estimate stays valid under assumptions ( |1.9| to ( |1.11[ ). We then write 



dA^OIlogel [1-C 



,log|logg| 
|loge| 



+d^ in (a J )F(a ] ) > C\dj\ g z { aj ) |loge| 1-C 



,log|logg| 
|loge| 



(A.19) 



-\F(a 3 )\ 



On the other hand, we have from the definition (3.30) 

IKI-i^l < Ce\loge\n{ /2 < e|loge| 



for any j S J out . The result (3.35) follows using (A.19). 

We turn to the energetic cost of the vortices in J ln . First we deduce from (A.16) that 

TF / — ' 1 



F(a j )=F lt (a j ) + 



v e| loge| 

for any j £ J in . It is then straightforward (recall that Hi ^> log | log £ 1 1 log el" 1 ) from the computation 
(A. 5 1 to obtain 



F( aj ) < - 



C 

e\ logej 1 /" 



< 
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for any j € J in such that \a-j\ < 1 — e\ logel 1 / 2 . On the other hand, if |oj| > 1 — e| logel 1 / 2 , g 2 (dj) > 
C(£|loge|) _1 and thus, using (A. 19) 

\F( aj )\ «<? 2 ( aj ). 

We deduce that for any j e J_ 



\g*{ aj ) |loge| + djt in (a j )F(a j ) > C\d 3 \g 2 (a 3 )\\oge\ > C\d^ 



1/2 



using (A. 15). On the other hand, if j 6 J + 
1 



2 5 (%) |loge| |di| +d^ ia (a j )F(a j ) = |d,-| I -<? >j) |loge| + &,(a i ).F(a i ) 



and thus (3.37) follows from Propositions A.l and A. 2 in the case where j € J+. We also use £i n < 1 



There remains to show that (3.36) holds when j 6 J*, which is a consequence of (A. 13) and (A. 18) once 
one has recalled that £i n (a.j) = 1 for any j e J,. □ 



Appendix B : Useful results from fCRYj 



In this appendix we gather for the convenience of the reader several results of jCRY] that are essential 
ingredients of the proof of our main results. 

We begin with estimates of the GP minimizer. The following is Proposition 2.2 of [CRYj . 

Proposition B.l (Exponential smallness of ^ GP inside the hole). 

As e — > and for any f € B, 



I^Mf^C^lloger 1 expj-^-^j. 

Moreover there exists a strictly positive constant c such that for any C(e 7 / 6 ) < r < i?h — 0(e 7 ^ 6 ) : 

|* GP (f)| 2 < epilogs!" 1 exp{--^}. 

We also employ several useful properties of g, starting with Proposition 2.4 of [CRY . 

Proposition B.2 (Preliminary estimates for g). 

As e — > and for any ui € Z such that \oj\ < 0(e^ 1 ), 



\9 2 -P TF \\ LHA) -0(l), 



< \\P TFI 



\\9\\l~(A) S \\P" \\ L co {A) il + C(Ve|loge| 
We next state the exponential decay of g |CRY[ Proposition 2.5]. 



Proposition B.3 (Exponential smallness of g inside the hole). 

As e — > and for any r G A 



g 2 (r) < Ce -1 |log£| -1 exp 



1-r 2 



Moreover there exists a strictly positive constant c such that for any r < R\ L — 0(e 7 / 6 ) ; 



>( r )<C £ -i|log £ |-iexp{-^} 



(B.l) 



(B.2) 



(B.3) 



(B.4) 



(B.5) 
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Finally, it is very useful to know that g 2 is very close to p TF in a L°° sense (CRY1 Proposition 2.6]. 

Proposition B.4 (Pointwise estimate for g). 

As e — > and for any uj £ Z with \w\ < 0(e^ 1 ) 

\g 2 (r) - p TF (r)\ < Ce 2 \ loge| 2 (r 2 - R 2 )- 3 / 2 p TF (r) (B.6) 
for any re V TF such that r > R h + 0(e 3/2 \ loge| 2 ). 
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